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ABSTRACT 


The  demand  for  pathways  to  fabricate  and  methodologies  to  understand  3D  structures 
with  advanced  functionalities  has  increased  significantly  over  the  last  years.  Materials  exhibiting 
three-dimensional  structure  with  characteristic  length  scales  ranging  from  nanometers  to 
micrometers  have  extraordinary  potential  for  emerging  optical  and  thermal  applications  and  are 
of  great  current  interest  in  providing  new  functionalities  for  a  host  of  applications. 

The  thesis  is  divided  into  two  primary  parts.  The  majority  of  the  thesis  is  devoted  to 
understanding  light-matter  interaction  in  a  new  class  of  devices  termed  Resonant  Guided  Wave 
Networks  (RGWNs).  We  describe  how  the  geometrical  properties  of  the  network  of  waveguides 
are  tuned  in  conjunction  with  the  materials  properties  to  realize  a  collective  optical  response, 
contrary  to  other  dispersive  photonic  materials  where  the  response  is  limited  to  its  discrete 
elements.  In  particular,  we  demonstrate  a  simple  and  efficient  fabrication  of  MIM  based  complex 
3D  structure  that  offers  a  novel  approach  to  optical  dispersion  control  based  on  resonant 
structures.  This  is  followed  by  a  description  of  experiments  that  probe  the  optical  response  of  the 
structures.  Of  particular  interest  are  surface  electromagnetic  modes  known  as  surface  plasmon- 
polaritons  (SPPs). 

In  the  second  part  of  this  thesis,  we  report  measurements  and  modeling  of  thermal 
conductivity  in  periodic  three-dimensional  dielectric  nanostructures:  silicon  inverse  opals.  Such 
structures  represent  a  three-dimensional  “phononic  crystal”  but  affect  heat  flow  instead  of 
acoustics.  The  thermal  conductivity  of  inverse  opal  films  are  relatively  low,  -0.6-1. 4  W/mK  at 
300  K  which  is  due  to  macroscopic  bending  of  heat  flow  lines  in  the  structure.  The 
corresponding  material  thermal  conductivity  is  in  the  range  5-12  W/mK  and  has  an  anomalous  - 
T18  dependence  at  low  temperatures,  distinct  from  the  typical  ~T3  behavior  of  bulk 
polycrystalline  silicon.  Using  phonon  scattering  theory,  we  show  such  dependence  arising  from 
coherent  phonon  reflections  in  the  inter-grain  region.  This  is  consistent  with  an  unconfirmed 
theory  proposed  in  1955.  The  low  thermal  conductivity  is  significant  for  applications  in 
photonics  where  they  imply  the  possibility  of  significant  temperature  rise  even  at  relatively  low 
optical  absorption  and  in  thermoelectrics,  where  they  suggest  the  possibility  of  enhancement  in 
the  figure  of  merit  for  polysilicon-based  devices. 
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PARTI 

PLASMONS  IN  COMPLEX  3D-STRU CTURE  (RGWNs)* 

*This  work  was  done  in  collaboration  with  Stanley  Burgos  from  Atwater  Research 
group  at  Caltech  and  Bok  Y.  Ahn  from  Lewis  Research  Group  at  School  of 
Engineering  and  Applied  Sciences,  Harvard. 
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CHAPTER  1 


INTRODUCTION  TO  SURFACE  PLASMONS 

1.1  Light-matter  interaction 

Light-matter  interactions  have  fascinated  mankind  for  thousands  of  years,  and  are  the 
basis  of  some  of  the  most  famous  scientific  discoveries  to  date.  With  the  continued 
miniaturization  and  integration  of  photonic  circuits,  recent  advances  are  becoming  greatly 
limited  by  diffraction.  In  classical  optics,  diffraction  restricts  the  localization  of  electromagnetic 
waves  into  nanoscale  regions  much  smaller  than  the  wavelength  of  light  in  the  material1. 
However,  the  electromagnetic  response  of  metals  has  led  to  the  development  of  the  emerging  and 
fast  growing  research  field  of  plasmonics.  Plasmonics  opens  a  path  for  controlling  light-matter 
interactions  on  the  subwavelength  scale,  serving  as  a  bridge  between  diffraction-limited  optics 
and  nanoscale  devices2. 

The  strong  interaction  between  light  and  metal  has  been  exploited  for  thousands  of  years, 
most  recognizably  in  the  art  of  creating  brightly  colored  stained-glass  window  panels  (Figure 
1.1a)  by  introducing  metallic  salts  in  otherwise  transparent  glass  matrix.  One  of  the  famous 
example  is  the  “Lycurgus  Cup”3  dating  back  to  the  4th  century  Rome  (Figure  1.1b).  Viewed 
normally  by  reflected  light  the  surface  is  opaque  olive  green,  but  when  backlit  the  cup  is  a 
brilliant  translucent  red.  This  fascinating  observation  is  due  to  the  inclusion  of  gold  and  silver 
particles  within  the  glass  matrix,  which  strongly  scatter  green  light  but  transmit  red. 
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a 
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Figure  1.1  -  (a)  Stained  glass  window  in  Rome  with  its  brilliant  colors,  (b)  The  Lycurgus 
cup  from  4th  century  Rome.  The  glass  consists  of  silver  and  gold  nanoparticles.  Image 
reproduced  from  Reference  3. 


Fundamental  research  and  development  of  metal  based  structures  and  devices  have 
attracted  a  number  of  groups  in  recent  years  to  manipulate  light  at  length  scales  smaller  than  the 
wavelength  of  light  by  employing  surface  plasmons  at  metal-air  and  metal-dielectric  interfaces. 
In  this  chapter,  I  review  the  electromagnetic  properties  of  metals  in  the  optical  frequency  regime. 
An  exploration  into  light-metal  interactions  will  then  be  presented,  building  up  to  studying  the 
properties  of  surface  plasmons  in  order  to  understand  how  they  can  be  employed  to  control  and 
manipulate  light  at  the  metal-dielectric  interface. 


1.2  Optical  properties  of  metals 

Metals  constitute  majority  of  the  periodic  table  and  have  been  studied  for  thousands  of 
years.  In  the  most  simple  description,  metals  can  be  thought  of  as  a  sea  of  free  electrons 
oscillating  within  a  lattice  of  fixed  ion  cores.  These  electrons  are  free  to  move  throughout  the 
metal  and  respond  to  any  applied  external  field.  However,  the  electron’s  response  to  an  optical 
field  is  entirely  dependent  on  the  frequency  response  of  the  metal’s  dielectric  constant4. 
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By  using  the  Drude  model  for  metals,  one  can  obtain  the  dielectric  constant  of  metals  as  a 


function  of  frequency  in  the  following  form 


eO)  = 


(1) 


where  u)p  is  the  plasmon  frequency  of  the  corresponding  bulk  metal,  and  r  is  the  electron 
relaxation  time  in  that  metal.  When  one  neglects  the  collision  (lossless  i.e.  co  »  r),  we  get 


6(oj)  = 


(2) 


Therefore,  metals  have  negative  dielectric  constants  at  frequencies  smaller  than  plasma 
frequency. 

At  low  (microwave  and  far-infrared)  frequencies,  many  metals  behave  like  a  perfect 
electrical  conductor  with  high  reflectivities  and  negligible  electromagnetic  field  penetration. 
Hence,  metals  are  not  the  first  materials  that  come  to  mind  for  designing  optical  elements. 
However,  above  certain  frequency,  the  driving  field  strongly  couples  into  the  longitudinal 
oscillations  of  the  bulk  electrons  within  the  metal4.  This  is  known  as  the  plasma  frequency  of  the 
metal  and  is  given  by 


where  n  is  the  density  of  electrons  within  the  metal  and  m*  is  the  effective  mass  of  the  electrons. 

For  noble  metals  of  interest  like  gold,  silver,  and  copper,  the  optical  frequency  regime 
lies  below  their  plasma  frequency  giving  them  a  negative  dielectric  constant  as  can  been  seen 
from  Equation  2.  When  such  metals  meet  another  material  with  a  positive  permittivity,  an 
electromagnetic  wave  can  be  bounded  by  their  interface  according  to  the  Maxwell  equations5. 
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1.3  Surface  Plasmon  Polaritons  (SPPs) 


As  discussed  in  the  previous  section,  the  negative  dielectric  constant  of  metals  at  optical 
frequencies  leads  to  existence  of  longitudinal  charge-density  fluctuations  along  a  metal-dielectric 
interface.  These  coherent  oscillations  of  free  electrons  are  referred  to  as  surface  plasmons  (SPs)6. 
SPs  can  couple  with  electromagnetic  waves  resulting  in  a  freely  propagating  surface  plasmon 
polaritons  (SPPs)  along  metal  surfaces  to  localized  surface  plasmons  on  metallic  nanoparticles 
(Figure  1.2).  Propagating  SPPs  allows  significant  reduction  in  effective  wavelength  and  a 
corresponding  significant  increase  in  spatial  confinement  and  local  field  intensity,  and  will  be  the 
focus  of  this  thesis. 


Dielec  trie 


Figure  1.2  -  Schematic  diagrams  illustrating  the  charge  and  electric  field  associated  with  (a)  a 
propagating  surface  plasmon  on  an  extended  metal  film  and  (b)  a  localized  surface  plasmon  on  a 
metal  nano-particle.  Image  reproduced  from  Reference  7. 

1.3.1  Literature  review 

The  first  scientific  observation  of  surface  plasmons  date  back  to  1902  by  Robert  W. 
Wood9.  In  1904,  J.C.M.  Garnett  first  describes  the  bright  colors  observed  in  metal  doped 
glasses10  employing  the  contemporary  Drude  model  of  free-electron  metals  and  the 
electromagnetic  properties  of  small  spheres  as  derived  by  Lord  Rayleigh11.  Shortly  thereafter,  in 
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1908,  Gustav  Mie  presented  a  general  formulation  for  the  scattering  of  light  from  spherical 
surfaces,  including  the  particular  case  of  colloidal  nanoparticles  of  varying  size  .  The  theoretical 
description  of  surface  plasmons  on  thin  films  was  first  described  in  terms  of  electron  energy  loss 

i  o 

spectroscopy  measurements  in  1957,  by  Ritchie  .  In  1968,  nearly  seventy  years  after  Wood’s 
original  observations,  Ritchie  and  coworkers  describe  the  anomalous  behavior  of  metal  gratings 
in  terms  of  surface  plasmon  resonances  excited  on  the  gratings14. 

A  major  experimental  milestone  in  the  study  of  surface  plasmons  was  made  in  1968  when 
Andreas  Otto15  as  well  as  Erich  Kretschmann  and  Heinz  Raether16  developed  methods  for  the 
optical  excitation  of  surface  plasmons  on  metal  films,  making  experiments  on  surface  plasmons 
easily  accessible  to  many  researchers.  As  the  field  continued  to  develop,  and  the  importance  of 
the  coupling  between  the  oscillating  electrons  and  the  electromagnetic  field  become  more 
apparent,  Stephen  Cunningham  and  his  colleagues  coined  the  term  “surface  plasmon-  polariton” 
(SPP)  in  1974  .  All  these  findings  have  set  the  stage  for  the  current  surge  in  surface  plasmon 
research  in  recent  years. 

1.3.2  Properties  ofSPPs 

The  electric  field  associated  with  SPPs  is  strongly  confined  (sub-wavelength)  and  decays 
exponentially  with  increasing  distance  from  the  interface  (Figure  1.3a)4.  Many  of  the  appealing 
features  of  SPPs  come  from  this  confinement:  The  electric  field  strength  in  the  vicinity  of  the 
interface  is  enhanced,  providing  possibilities  to  enable  a  wide  range  of  practical  applications 
including  light  guiding  and  manipulation  at  nanoscale  '  ,  biodetection  at  single  molecule 

23  25  26  27 

level  ,  enhanced  optical  transmission  through  subwavelength  apertures'  ’  ,  high  resolution 

98  30  31  39 

optical  imaging  below  the  diffraction  limit'  '  ,  and  solar  cells  ’  .  Recently,  considerable  focus 
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has  been  directed  to  exploit  the  dispersion  characteristics  of  SPPs  traveling  in  confined  metallo- 
dielectric  spaces  ’  to  create  photonic  materials  with  artificially  tailored  bulk  optical  properties, 
which  will  be  the  focus  of  this  thesis. 


Figure  1.3  -  (a)  Electric  field  distribution  perpendicular  to  the  surface.  The  SPP 
electromagnetic  field  decays  exponentially  with  increasing  distance  into  each  medium  from 
the  interface,  (b)  Dispersion  relation  for  a  surface  plasmon,  light  line  represents  dispersion 
for  light  in  a  vacuum. 


The  short-wavelength  SPPs  enable  the  fabrication  of  nanoscale  optical  integrated  circuits, 
in  which  light  can  be  guided,  split,  filtered,  and  even  amplified  using  plasmonic  integrated 
circuits  that  are  smaller  than  the  optical  wavelength.  Since  photons  can  be  coupled  to  plasmons 
and  vice  versa  and  the  frequency  range  of  SPPs  extends  over  the  infrared  and  visible  spectrum, 
they  also  provide  a  new  platform  for  optical  communication  .  However,  their  decay  due  to  both 
scattering  from  the  surface  imperfections  and  ohmic  losses  within  the  metal  has  limited  their 
application.  With  the  advances  in  micro-nano  fabrication  technologies,  there  has  been  significant 
improvement  in  SPP  properties  particularly  the  choice  of  metal  (dispersion),  metal  thickness,  and 


excitation  frequency36. 
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1.3.3  Dispersion  relation 


The  electromagnetic  field  associated  with  SPP  at  a  dielectric/metal  interface  can  be 
obtained  by  solving  Maxwell’s  equations4.  From  the  solution  one  obtains  SPP’s  dispersion 
relation,  that  is,  the  SPP  energy  (or  frequency)  as  a  function  of  its  wave- vector  (Figure  1.3b), 


edem 
ed  +  em 


=  kx  +  i  k 


X 


(4) 


where  ed  and  em  are  the  frequency  dependent  dielectric  function  of  the  dielectric  and  metal 
respectively. 

The  SPP  dispersion  relation,  co  (k)  plotted  in  Figure  1.3b,  shows  that  for  small  wave 
vectors  the  plasmon  is  close  to  that  of  the  light  line,  but  always  to  the  right  hand  side.  Therefore, 
the  plasmon  has  similar  properties  to  the  optical  field  but  always  has  a  greater  momentum  and  is 
non-radiative.  This  increased  momentum  is  associated  with  the  binding  of  the  mode  to  the 
surface.  For  larger  wave  vectors  the  surface  plasmon  dispersion  tends  towards  a  maximum  value 
of  c Up /  s/2  ,  and  resembles  that  of  the  bulk  plasmon  and  is  radiative,  the  energy  becoming 
independent  of  wave  vector. 

Since  SPPs  do  not  directly  couple  to  an  optical  field  they  must  be  created  via  some 
intermediate  step1.  This  requires  either  converting  the  incident  light  into  an  evanescent  mode  by 
means  of  total  internal  reflection,  or  by  scattering  the  light  to  higher  momentum  via  a  grating  or 
scattering  site.  Basically,  some  mechanism  must  be  introduced  to  effectively  increase  the  wave 
vector  of  the  incident  light  in  order  to  achieve  momentum  matching  between  the  surface 
plasmons  and  the  free-space  photons  at  a  given  energy  as  evident  from  Figure  1.3b.  While 
discussing  the  dispersion  relation  it  is  of  worth  to  discuss  the  spatial  extent  of  the  fields 
associated  with  the  SPP  mode. 
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1.3.4  Propagation  length  and  skin  depth 

Due  to  the  imaginary  part  of  em,  the  propagation  distance  of  SPPs  is  limited  due  to 
absorption  in  the  metal.  The  finite  propagation  length  is  accounted  for  by  the  imaginary  part  of 
k$pP.  As  the  SPP  mode  propagates,  the  electric  field  amplitude  decays  exponentially  as 
exp(—kxx)  and  therefore  the  intensity  as  exp(—2kxx).  We  can  define  the  propagation  distance 
L  as  the  1/e  decay  length  of  the  field  intensity  along  x, 

L=^77«-(^ (5) 

2 kx  G)  V  £d^m  '  £ m 

Here  e'm  and  are  the  real  and  imaginary  parts  of  the  dielectric  function  of  the  metal  such  that 
e  —  e'  +  ie" 

The  propagation  distance  depends  on  the  dielectric  constant  of  the  metal  and  the 
operating  frequency  ranging  from  10pm  in  the  visible  to  millimeter  in  the  infrared  for  metals 
with  low  adsorption  loss  such  as  silver,  gold4.  The  field  penetration  (plasmon  skin  depth)  within 
the  dielectric  and  the  metal  are  of  the  order  of  tens  to  hundreds  of  nm  as  shown  in  Figure  1.3a. 

1.4  Metal-Insulator-Metal  (MIM)  waveguide 

As  mentioned  earlier,  surface  plasmons  on  a  planar  metal-dielectric  interface  allow 
significant  reduction  in  the  effective  wavelength  and  a  corresponding  significant  increase  in 
spatial  confinement  and  local  field  intensity  of  the  light.  Placing  a  second  metal  layer  above  the 
dielectric,  one  can  achieve  a  metal-insulator-metal  (MIM)  waveguide  with  extremely  high  modal 
confinement  as  shown  in  Figure  1.4.  This  significantly  reduces  the  decay  of  optical  field  by 
confining  the  electromagnetic  field  in  the  dielectric  region. 
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Provided  the  thickness  of  the  metal  film  exceeds  the  plasmon  skin  depth,  oscillations  at 
each  metal-dielectric  interface  are  decoupled,  and  independent  surface  plasmon  modes  at  each 
metal-dielectric  interface  are  sustained. 

a  b 


Figure  1.4  -  Schematic  showing  the  electric  field  lines  of  a  surface  plasmon  (a)  at  a 
metal-dielectric  interface,  and  (b)  for  a  metal-insulator-metal  waveguide.  Image 
reproduced  from  Reference  37. 

As  seen  in  Figure  1.4,  the  MIM  geometry  significantly  confines  the  electromagnetic 
fields  in  the  dielectric  region  more  than  a  single  metallodielectric  interface.  By  using  dielectrics 
with  higher  refractive  index,  one  can  decrease  the  element  size  of  the  plasmonic  device  to  the 
nanometer  scale  by  localizing  light  within  thinner  insulator  thickness.  In  addition,  the  output  of 
the  device  can  be  tuned  across  the  visible  and  infrared  region  of  the  electromagnetic  spectrum 
by  replacing  the  insulator  layer  with  an  active  material.  By  changing  each  of  the  materials  in 
the  three  layer  stack  (Figure  1.4b),  as  well  as  their  relative  thicknesses,  the  dispersion 
properties  of  the  device  as  well  as  the  resonance  and  filtering  properties  can  be  tuned  over  a 
wide  range  of  the  electromagnetic  spectrum.  In  addition,  by  changing  the  dimensions  of  the 
structures,  as  well  as  the  distance  between  them,  we  can  selectively  access  a  large  range  of  k- 
vectors  or  wavelengths  that  are  present  within  the  structure. 


dielectric 


metal 
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CHAPTER  2 


RESONANT  GUIDED  WAVE  NETWORKS  (RGWNs) 

2.1  Dispersive  photonic  materials 

Nanofabricated  systems  that  exploit  SPPs  demonstrate  potential  for  controlling  light- 
matter  interaction  on  subwavelength  scale,  enabling  nanophotonic  technologies  that  are 
unattainable  with  conventional  diffraction-limited  optical  systems^.  SPPs  are  able  to  channel 
light  efficiently  into  nanometer  scale  volumes,  leading  to  direct  modification  of  resonate 
frequency  dispersion  properties  (substantially  shrinking  the  wavelength  and  speed  of  light  for 
example),  as  well  as  field  enhancements  suitable  for  enabling  strong  nonlinear  effects4.  In  this 
thesis,  we  introduce  a  novel  approach  to  optical  dispersion  control  based  on  three-dimensional 
(3D)  metallo-dielectric  periodic  structures  by  exploiting  the  dispersion  characteristics  of  SPPs. 
Materials  exhibiting  three-dimensional  structure  with  characteristic  lengths  ranging  from 
nanometers  to  micrometers  have  extraordinary  potential  for  emerging  optical  applications  based 
on  the  regulation  of  light-matter  interactions  via  the  local  interaction  of  periodic  elements.  The 
majority  of  this  chapter  is  focused  on  the  design,  fabrication,  and  experimental  characterization 
of  these  resonant  based  plasmonic  structures. 

Recent  years  have  seen  a  rapid  growth  of  research  into  various  artificial  photonic  material 
concepts  that  have  defined  new  directions  to  control  optical  dispersion,  thus  exhibiting  superior 
electromagnetic  properties.  Most  of  the  work  has  been  dedicated  to  metamaterials  due  to  their 
new  and  seemingly  exotic  electromagnetic  response.  Metamaterials  are  artificially  engineered 
metallo-dielectric  structures  patterned  to  sizes  much  smaller  than  the  operating  wavelength,  and 
their  special  geometrical  arrangement  affects  the  permittivity  and  permeability  of  the 
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electromagnetic  field  to  achieve  light  propagation  and  manipulation  different  than  their 
constituent  materials,  in  a  way  not  achievable  through  natural  medium  ’  ’  .  These  engineered 
metamaterials  enable  a  host  of  novel  effects  such  as  negative  index  of  refraction,  optical 
cloaking,  and  sub -diffraction  imaging.  They  commonly  require  complex  multidimensional 
building  blocks,  which  are  difficult,  if  not  impossible,  to  achieve  via  standard  technology.  Also, 
use  of  metals  at  higher  frequency  can  be  problematic  due  to  their  losses.  As  an  alternative, 
researchers  investigated  into  periodic  structures  typically  composed  of  low  loss  dielectrics 
known  as  photonic  crystals  (PCs).  PCs  offer  unprecedented  control  of  light  propagation  in  a  way 
similar  to  how  the  periodicity  of  crystals  effect  the  motion  of  electrons  by  defining  allowed  and 
forbidden  energy  bands40'41.  PCs  are  the  optical  analogs  to  semiconductors  that  have  gained  a  lot 
of  interest  especially  due  to  their  ability  to  suppress  light  propagation  in  all  three  dimensions 
(complete  photonic  band  gap)  .  In  other  words,  PCs  provide  a  platform  to  engineer  light-matter 
interaction  at  will.  Another  novel  approach  to  light  manipulation  utilizes  high  Q  optical 
resonators  that  are  kept  in  close  proximity  to  one  another.  The  resonators  are  based  on  whispery 
gallery  modes  and  are  able  to  manipulate  light  propagation  through  evanescent  coupling  of  the 
electromagnetic  field43. 

In  this  chapter,  I  will  focus  on  a  fairly  new  approach  for  optical  dispersion  control  based 
on  resonant  guided  network  of  MIM  waveguides.  These  structures  were  recently  proposed  and 
are  referred  to  as  resonant  guided  wave  networks  (RGWNs)44. 
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2.2  Resonant  Guided  Wave  Networks  (RGWNs) 


RGWNs  offer  a  novel  approach  to  optical  dispersion  control  based  on  localized 
resonance  within  these  structures44.  RGWNs  are  normally  intersecting  MIM  waveguides  that 
serve  as  wave  splitting  elements  arranged  in  two-  or  three-dimensional  networks  (Figure  2.1a). 
The  resulting  network  resonances  give  rise  to  wave  dispersion  that  is  sensitive  to  the  network 
layout,  contrary  to  other  dispersive  photonic  systems.  RWGNs  have  demonstrated  potential  for 
highly  efficient  resonators,  band-gap  materials44  and  devices  with  programmable  optical 
transmission  properties45. 

We  will  first  introduce  the  properties  of  RGWN  building  block  and  investigate  the  wave 
dynamics  within  these  structures  using  the  finite-difference  time-domain  (FDTD)  method  with  a 
short  pulse  excitation  at  a  central  wavelength  of  1.5  pm.  A  resonant  guided  wave  network  is 
comprised  of  power  splitting  elements  connected  by  isolated  waveguides  (Figure  2.1a). 


2D  rgwn  3D  rgwn  Equal  Mode  Splitting  Plasmonic  Resonators 


Figure  2.1  -  (a)  Schematic  of  RGWN  building  block,  (b)  Equal  mode  splitting  in  RGWNs 
(left);  RGWNs  as  optical  resonators  (right).  Image  reproduced  from  Reference  44. 


The  power  splitting  elements  distribute  the  incoming  pulse  to  the  waveguides  connected 
to  the  junction  with  certain  phase  and  amplitude.  It  was  recently  found  that  a  cross  junction  of 
two  normally  intersecting  MIM  waveguides  exhibit  equal  splitting  of  optical  power  at 
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subwavelength  gap  sizes46.  The  subwavelength  MIM  junction  serves  as  a  nano-aperture 
(scattering  site)  for  the  incoming  pulse,  resulting  in  a  broad  spectrum  of  plane  waves  with 
significant  fraction  of  the  optical  power  coupling  into  the  waveguides  perpendicular  to  the  input 
waveguide.  This  strong  coupling  of  the  incoming  pulse  to  all  neighboring  waveguides  gives  this 
structure  an  optical  response  different  from  purely  dielectric  waveguides,  where  most  of  the 
power  is  transmitted  in  the  forward  direction.  The  nano-aperture  effect  is  limited  for  dielectric 
waveguides  by  their  half-wavelength  modal  cross-sections. 

In  addition  to  the  equal-power  splitting,  the  transmitted  and  the  reflected  pulses  are  out- 
of-phase.  The  pulses  undergo  interference  inside  the  network,  and  interfere  destructively  or 
constructively  depending  on  the  waveguide  length  and  network  topology.  Hence,  if  one 
appropriately  designs  the  junction  and  waveguide  components  as  well  as  the  network  layout  with 
subwavelegth  gap  size  and  intemode  distance  on  the  order  of  the  effective  wavelength,  a  small  2 
x  2  resonant-guided  wave  networks  can  function  as  an  ultracompact  optical  resonator  (Figure 
2.1b).  By  controlling  the  properties  of  each  splitting  element  and  every  waveguide,  the  network 
interference  and  resonances  can  be  determined  and  thus  also  the  optical  function  of  the 
network45.  Also,  the  local  network  topology  can  be  adjusted  to  operate  the  structure  over  broad 
range  of  frequencies. 

Now  if  we  consider  an  infinite  2D  RGWN  structure,  photonic  bandgaps  are  observed  due 
to  the  interference  effect  as  shown  in  Figure  2.2a.  In  other  words,  certain  wavelength  of  light 
would  enter  a  resonant  state  and  not  get  transmitted,  giving  this  structure  a  photonic  bandgap. 
The  photonic  bandgap  can  be  engineered  by  controlling  the  dispersion  properties  of  the  network: 
the  phase- shift  upon  splitting  which  is  controlled  by  the  width  of  the  MIM  gaps  and  the 
interference  which  is  dependent  on  the  distance  between  the  nodes. 
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Figure  2.2  -  (a)  Optical  density  of  states  (right)  for  an  infinite  2D  RGWNs  (left),  (b)  Optical 
density  of  an  infinite  3D  network.  Image  reproduced  from  Reference  44. 


The  wave  dynamics  is  the  same  for  3D  RGWNs,  however,  the  four-arm  cross  junction 
element  of  the  2D  network  is  replaced  by  a  six  arm  cross  junction  in  a  3D  network  (Figure  2.1a). 
It  has  been  demonstrated  that  six-way  equal  power  splitting  occurs  in  the  3D  coaxial  MIM 
waveguide  networks44.  Likewise  the  dispersion  properties  are  governed  by  the  local  network 
topology  giving  the  structure  tunable  photonic  bandgap  (Figure  2.2b). 


2.3  Fabrication  process 

2.3.1  Direct-ink  writing 

The  ability  to  pattern  materials  in  three  dimensions  is  critical  for  a  number  of  emerging 
technologies,  including  photonics47’48,  microfluidics49,  microelectromechanical  systems50, 
biomaterials51’52,  batteries53,  and  optoelectronics54.  I  specifically  have  focused  on  direct- writing 
fabrication  techniques,55  which  allows  one  to  design  and  rapidly  fabricate  complex  three- 
dimensional  structures  without  the  need  for  expensive  and  complicated  tools,  dies,  or 
lithographic  routes.  In  this  method,  a  concentrated  ink  is  extruded  through  a  tapered  cylindrical 
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nozzle  that  is  translated  onto  the  substrate  using  a  computer-controlled  3-axis  micropositioning 
stage  (Figure  2.3a).  Due  to  the  viscoelastic  ink  characteristics,  direct-ink  writing  enables  self- 
supporting  flexible,  stretchable  and  spanning  features  from  0.1pm-250pm55.  This  low-cost, 
large-area  fabrication  method  with  higher  flexibility  in  material  selection  overcomes  many 
challenges  that  exist  with  the  pre-existing  printing  techniques  for  fabrication  of  complex  3D 
structure. 

We  use  a  direct-written  3D  woodpile  structure  as  a  scaffold  to  fabricate  the  coaxial  3D 
RGWNs  (Figure  2.3b).  The  scaffold  is  composed  of  highly  concentrated  silver  nanoparticles  ink 
with  a  solid  loading  of  ~80%56.  The  scaffold  network  is  coated  with  thin  Ag  film  (-150  nm)  via 
electrodeposition,  and  then  coated  with  silicon  (-200  nm)  using  Si-CVD  followed  by  a  second 
Ag  eletrodeposition.  The  thicknesses  were  appropriately  chosen  to  ensure  coherent  plasmonic 
mode  propagation  and  the  resonant  behavior  using  FDTD  simulations.  The  schematic  of  the 
structure  is  shown  below  (Figure  2.3c-d). 
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Figure  2.3  -  (a)  Schematic  illustration  of  direct-write  assembly.  Image  reproduced 
from  Reference  56.  (b)  SEM  image  of  the  direct- written  3D  woodpile  scaffold.  Inset 
shows  the  magnified  view,  (c)  Schematic  of  the  unit  structure  for  3D  RGWN 
structure,  (d)  Schematic  of  the  cross-section  of  the  3D  RGWNs. 
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The  3D  scaffold  is  composed  of  silver  nanoparticles  along  with  some  polymer  additive  to 
promote  uniform  printing  contrary  to  pure  silver.  This  composite,  particle-based  structure  is  not 
ideal  for  plasmonic  properties.  Hence,  we  coat  the  structure  with  a  thin  conformal  silver  film. 
However,  when  printed,  the  silver  nanoparticles  are  not  yet  bonded  together,  and  thus  thermal 
annealing,  which  allows  removal  of  the  organic  additives  to  some  extent  and  promote 
nanoparticle  sintering,  is  performed  before  coating  with  a  thin  layer  of  conformal  silver.  The 
thermal  treatment  increases  the  electrical  conductivity  of  the  printed  3D  scaffold56.  This  is 
critical  for  the  uniform  and  conformal  electro-deposition  of  silver  film.  The  conductivity  of  the 
ink  as  a  function  of  annealing  temperature  is  shown  in  Figure  2.4. 


Figure  2.4  -  Resistivity  of  direct- written  silver  ink  films  as  a  function  of  annealing 
temperature.  Right  shows  top  view  SEM  images  of  the  silver  coating  on  annealed 
silver  ink  films. 

For  annealing  temperatures  below  250°C,  the  conductivity  of  the  silver  ink  is  not  high 
enough  for  silver  deposition.  However,  at  higher  annealing  temperatures  the  deposited  silver  gets 
rougher.  As  discussed  in  Section  1.3,  surface  plasmons  are  very  sensitive  to  surface 
inhomogeneties  as  SPPs  exist  and  propagate  very  close  to  the  interface.  Hence,  optimizing  the 
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surface  roughness  and  other  homogeneities  is  very  critical  while  fabricating  these  structures  as 
any  inhomogeneities  lead  to  scattering  losses  and  limited  propagation  in  addition  to  the  intrinsic 
absorption  loss  from  the  metal.  Also,  the  surface  morphology  of  the  electrodeposited  silver 
determines  the  surface  morphology  of  both  the  metal-dielectric  interface,  making  surface 
roughness  of  electrodeposited  silver  particularly  detrimental  to  device  performance. 

2.3.2  Thermal  stability  of  the  direct-written  structure 

Thermal  annealing  increases  the  mechanical  stability  of  the  structure  which  is  critical 
during  further  processing  at  elevated  temperatures.  However,  the  silver  scaffold  suffers  serious 
structural  degradation  due  to  grain  growth  and  voids  at  high  temperature  and  would  often  distort 
and  delaminate  from  the  substrate. 


Annealed  in  air 


Annealed  in  forming  gas 


Figure  2.5  -  Top  view  SEM  images  of  silver  ink  films  after  annealing  at  different 
temperature.  Inset  shows  the  cross-section  view.  Top  samples  were  annealed  in  air  and 
bottom  samples  in  forming  gas. 
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We  were  able  to  achieve  the  structural  integrity  even  at  higher  temperature  by  annealing 
the  structure  under  forming  gas  (5%H2  in  Ar)  as  opposed  to  air56  (Figure  2.5).  The  ideal 
annealing  condition  (250°C  for  3  hrs,  ramping  rate  2°C/min'1)  was  determined  experimentally  for 
optimized  surface  roughness  of  the  electroplated  silver  (Section  2.3.4 )  on  the  printed  3D 
scaffold. 

2.3.3  Choice  of  metal 

As  discussed  earlier,  metallic  components  are  invariably  accompanied  by  losses  that  can 
be  detrimental  to  the  SPPs  properties  thus  limiting  the  feasibility  of  many  novel  plasmonic 
applications.  Hence,  the  choice  of  metal  is  very  critical  while  fabricating  these  plasmonic 
structures  to  achieve  better  device’s  performance.  Owing  to  their  relatively  lower  loss,  noble 
metals  such  as  gold,  silver,  and  copper  are  conventionally  used  as  metallic  components  at  optical 
frequencies. 

The  SPP  property  is  entirely  dependent  on  the  frequency  response  of  the  metal’s 
dielectric  constant.  The  real  part  of  the  dielectric  function  (e'm)  describes  the  strength  of  the 
surface  plasmons  (polarization)  induced  by  an  external  electric  field,  and  the  imaginary  part  {elf) 
describes  the  losses  encountered.  The  metal  of  choice  is  associated  with  large  negative  values  of 
e'm  and  small  values  of  e'f  .  Hence,  silver  is  the  metal  of  choice  for  plasmonics  in  the  near- 
infrared  region57. 

2.3.4  Electrodeposition  of  silver 

After  the  annealing  condition  was  optimized  @  250C  for  3  hrs,  the  3D  woodpile  scaffold 
was  coated  with  a  very  thin,  conformal  silver  film  via  electroplating.  Electroplating  is  a  process 
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in  which  metal  ions  in  a  solution  are  moved  by  an  electric  field  from  anode  to  coat  an  object 
connected  on  a  lower  electric  potential  known  as  cathode.  Silver  was  electrodeposited  using  a 
three-electrode  electrochemical  cell  consisting  of  a  platinum  counter  electrode  (anode),  an 
Ag/AgCl  reference  electrode,  and  sample  as  the  working  electrode  (cathode)  as  shown  in  the 
schematic  (Figure  2.6).  The  anode  and  cathode  are  both  connected  to  an  external  supply  of 
current.  The  anode  is  connected  to  the  positive  terminal  of  the  supply,  and  the  cathode  (sample  to 
be  plated)  is  connected  to  the  negative  terminal.  When  the  current  is  switched  on,  the  metal 
cations  in  the  electrolyte  move  to  the  cathode,  they  are  reduced  by  gaining  electrons  to  deposit 
the  article  with  their  metallic,  zero  valence  state. 


Figure  2.6  -  Schematic  of  the  3-electrode  electrodeposition  setup 

Commercial  plating  solution  of  silver  (Ag)  (Techni  Silver  CY-LESS  II  RTU)  purchased 
from  Technic,  Inc  was  used  as  the  electrolyte  and  patterned  gold  was  used  as  contact  to  the  3D 
silver  ink  scaffold  for  the  working  electrode.  The  sample  was  immersed  in  Isopropyl  alcohol 
(IPA)  prior  to  immersion  into  the  aqueous  plating  solution.  IPA  lowered  the  surface  tension 
allowing  full  penetration  of  the  electrolyte  into  the  complex  3D  structure.  The  electrochemical 
reaction  increased  in  efficiency  when  the  plating  solution  was  agitated  and  heated  to  40°C  during 
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electrodeposition.  At  40°C,  we  also  observed  reduction  in  the  grain  size  of  the  electrodeposited 
silver,  resulting  in  smoother  surface  morphology. 

After  electroplating,  we  observed  non-conformal  silver  deposition  for  constant  potential, 
more  deposition  in  areas  of  sample  closer  to  the  counter  electrode  due  to  diffusion  limitation.  To 
circumvent  the  diffusion  problem,  we  employed  pulsed  electrodeposition  technique  and 
optimized  it  experimentally.  The  deposition  potential  was  applied  for  a  short  interval  (-0.55V  for 
0.5seconds)  with  off-time  for  a  longer  interval  (0  V  for  3  seconds).  The  off-time  allowed  the 
metal  ions  in  solution  to  uniformly  diffuse  throughout  the  complex  3D  scaffold.  Using  this 
method,  we  observed  uniform  and  conformal  deposition  throughout  the  structure  (Figure  2.7b). 


Figure  2.7  -  SEM  image  of  silver  coating  on  ID  direct- written  line  templates  when 
employing  (a)  a  constant  potential  and  (b)  a  pulsed  potential 
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Figure  2.7  shows  SEM  images  after  electrodeposition  on  silver  lines.  The  dark  region 
shows  the  area  deposited  with  silver.  We  can  clearly  see  that  the  deposition  is  non-uniform  for 
constant  potential.  The  thickness  of  the  silver  film  was  controlled  by  the  number  of  iterations  of 
the  electrodeposition  sequence  (one  sequence  includes  one  deposition  potential  and  one  off- 
time). 

For  thick  layers  of  silver,  we  obtained  very  smooth  coatings  (Figure  2.7b).  However, 
increase  in  feature  size  of  the  structure  results  in  higher  order  modes  within  the  structure 

co 

resulting  in  modal  cutoff  in  the  visible  and  near-infrared  wavelength  range  .  Hence,  we  want 
very  thin  (~5  x  skin  depth)  and  smooth  silver  layer.  Smooth  surface  is  hard  to  achieve  for  a  thin 
layer  since  electrodeposition  is  conformal  and  follows  the  surface  morphology  of  the  underlying 
layer  as  shown  in  Figure  2.8a.  We  employed  electropolishing  to  reduce  the  surface  roughness. 

2.3.5  Electropolishing  of  silver 

By  far,  the  most  effective  way  to  electropolish  silver  is  the  “peroxide-cyanide  bomb”59’60. 
This  process  is  however  highly  volatile  due  to  the  presence  of  hydrogen  peroxide.  The  process 
should  only  be  done  in  a  dedicated  glove  box  designed  to  handle  toxic  gases  because  of  the 
explosive  nature  of  the  process  and  large  quantities  of  hydrogen  cyanide  that  is  produced. 
(Hydrogen  cyanide  is  very  toxic  and  should  always  be  handled  inside  a  glove  box).  Also,  the 
solution  should  be  made  fresh  every  time  to  achieve  good  polishing  results. 

In  this  work,  we  employed  a  much  safer  and  non-toxic  electrolytic  solution  for  anodically 
producing  a  smooth  silver  surface  through  electropolishing.  The  principle  carrier  electrolyte  was 
thiourea,  which  is  then  stabilized  using  larger  quantities  of  a  reducing  sugar  (glucose)61’62.  In 
addition,  the  electrolyte  uses  an  activating  acid  (methane  sulfonic  acid)  to  activate  the  thiourea 
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solution  and  urea  to  promote  uniform  dissolution  of  elements  and  provide  solution 
conductivity61.  Glycerin  is  used  as  a  surfactant  to  reduce  the  surface  tension  between  the  sample 
and  the  electrolyte  to  facilitate  uniform  electropolishing.  The  electropolishing  solution  is 
environmentally  safe,  operates  over  a  broad  range  of  current  density,  fast,  and  most  importantly 
is  functional  on  irregular  shapes  and  complex  designs. 

The  electrolyte  consists  of 

a.  thiourea:  100  grams  per  liter  of  solution; 

b.  glucose:  100  grams  per  liter  of  solution; 

c.  urea:  50  grams  per  liter  of  solution; 

d.  methane  sulfonic  acid:  20  milliliter  per  liter  of  solution;  and 

e.  glycerin:  50  milliliter  per  liter  of  solution. 

First,  a  volume  of  millipore  water  is  placed  in  a  vessel  equal  to  one-half  of  the  desired 
final  solution  volume.  The  water  is  warmed  to  a  temp  of  ~38°C  on  a  hot  plate  with  mechanical 
stirring.  Thiourea  and  glucose  is  then  added  to  the  water  slowly  which  cools  the  solution.  It  is 
desirable  to  maintain  the  temperature  of  ~38°C  while  adding  the  thiourea  and  glucose  for  better 
results.  The  solution  temperature  falls  to  ~30°C  after  thiourea  and  reducing  sugar  have  all 
dissolved.  Methane  sulfonic  acid  is  added  slowly  followed  by  glycerin.  The  sample  is  made 
anode  (working  electrode)  with  platinum  as  counter  electrode  as  shown  in  Figure  2.6.  We  do  not 
use  a  reference  electrode  here. 

The  electrolyte  temperature  is  maintained  at  ~30°C  with  constant  mechanical  agitation  of 
450  rotations  per  min  throughout  the  polishing.  We  employed  a  pulsed  potential  in  the  range  5V- 
8V  to  facilitate  uniform  electropolishing.  SEM  was  used  to  characterize  the  quality  of  the 
electropolishing.  Best  results  were  obtained  when  the  deposition  potential  was  applied  for  a  very 
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short  interval  (7V  for  40.0  ms)  and  off-time  for  a  lager  interval  (0  V  for  16.0  seconds)  with  black 


finish  observed  for  high  current  or  slow  agitation  (Figure  2.8b). 


Optimized  condition: 

■  Voltage  (Pulsed):  7.0  V 

■  On-time:  40.0  ms 

■  Off-time:  16.0  s 


Figure  2.8  -  Top  view  SEM  images  of  silver  coating  after  (a)  electroplating  of  very  thin 
silver  layer  (-130  nm)  (b)  electropolishing.  Inset  shows  the  magnified  view. 

After  electropolishing,  the  surface  is  still  not  quite  as  smooth  as  we  want  (Figure  2.8b 
Inset).  Hence,  we  further  optimized  the  surface  roughness  and  overall  quality  of  the  silver 
coating  via  thermal  annealing. 

2.3.6  Thermal  annealing  and  its  effect  on  the  silver  coating 

Post- thermal  annealing  is  a  popular  technique  used  in  plasmonics  to  improve  metal’s 
optical  and  electrical  quality64.  The  quality  of  the  silver  coating  was  significantly  improved  by 
thermal  annealing  at  200°C  in  forming  gas  (5%H2  in  Ar)  for  20  mins  (ramp  rate  2°C/min). 
Figure  2.9a  shows  a  significant  increase  in  the  reflectance  spectra  of  the  annealed  silver  films  in 
comparison  to  as  deposited  (electroplated  +  electropolished).  After  thermal  annealing,  the 
dimensions  of  the  silver  increase  significantly,  suggesting  sintering  and  coarsening  processes  are 
active.  This  welding  of  the  particles  and  disappearance  of  the  voids  and  discontinuities  (Figure 
2.9c)  responsible  for  scattering,  results  in  an  overall  film  quality  improvement.  The  average 
electrical  resistivity  of  the  annealed  silver  films  measured  by  four-point  probe  is  4.83  x  10"6  Q  ■ 
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cm.  This  value  is  an  order  of  magnitude  lower  than  before  annealing  (1.57  x  10'5  Q  ■  cm)  and 
approaches  near  bulk  silver  resistivity  (~  1.59  x  10  6  Q  •  cm). 


Figure  2.9  -  (a)  Reflectance  spectra  of  a  flat  silver  coating  after  annealing.  Top  view  SEM 
images  of  silver  coating  after  (b)  electroplating  +  electropolishing  and  (c)  thermal  annealing. 


Atomic-force  microscopy  was  used  to  quantify  the  surface  morphology  evolution  of  the 
deposited  silver  films  during  the  optimization.  All  AFM  images  were  obtained  with  an  Asylum 
Research  MFP-3D  system.  The  samples  were  scanned  using  tapping  mode,  typically  over  a  5x5 
pm'  area  at  a  scan  rate  of  1  Hz.  AFM  scans  show  a  significant  decrease  in  the  surface  roughness 
of  the  silver  films  (Figure  2.10). 

After  optimizing  the  silver  coating,  we  deposit  a  very  thin  layer  of  amorphous  silicon  on 
the  3D  silver  template  via  CVD. 
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Figure  2.10  -  AFM  profile  showing  the  reduced  roughness  of  flat  silver  films  after 
(a)  direct-ink  writing,  (b)  electroplating  +  electropolishing  and  (c)  thermal  annealing.  RMS 
roughness:  (a)  ~70  nm,  (b)  -15  nm  and  (c)  ~5  nm 


2.3. 7  Chemical  Vapor  Deposition  ( CVD)  of  silicon 

A  thin  conformal  film  of  high-quality,  low-roughness  amorphous  Si  was  deposited 
around  the  3D  template  using  a  static  chemical  vapor  deposition  (CVD)  system65'90  with  disilane 
(Si2H6,  98%,  Gelest)  as  a  silicon  source  at  350°C  (~50  mbarr,  heating  rate  8°C  min"1)  inside  a 
furnace.  A  static  CVD  system,  where  the  reacting  gas  was  confined  in  a  cell  together  with  the 
sample,  was  utilized  since  low  pressure  CVD  systems  show  negligible  growth  rates  at  this 
temperature66.  The  decomposition  time  controlled  the  thickness  of  the  silicon  film.  A  5  hr  run  at 
350°C  deposited  -70  +  5  nm  of  a-silicon  which  was  later  confirmed  using  SEM.  It  was  also 
observed  that  higher  pressures  did  not  significantly  increase  the  growth  rate,  while  samples 
prepared  at  higher  pressures  had  a  tendency  to  peel  off  the  substrate.  Because  of  the  static  CVD 
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setup,  it  was  necessary  to  grow  the  silicon  in  several  cycles  rather  than  one  longer  run,  because 
the  reaction  products  are  not  removed  from  the  chamber  and  eventually  the  reaction  reaches 
equilibrium.  For  silicon  thickness  of  ~200nm,  three  cycles  were  performed  and  was  later 
confirmed  using  SEM. 

2.3.8  Electroless  plating  of  silver 

The  amorphous  silicon  is  nonconductive.  There  are  various  methods  to  coat  silver  on  a 
nonconductive  substrate  such  as  sputtering,  plating,  chemical  vapor  deposition,  etc.  All  of  these 
techniques  have  their  own  pros  and  cons.  The  most  commonly  used  method  in  plasmonics  is 
plating,  due  to  its  ease  of  operation,  capability  to  coat  thick  and  isotropic  coating,  and  low  cost. 
In  this  work,  we  employ  electroless  plating  to  achieve  conformal  coating  of  the  third  metal  layer 
onto  the  complex  3D  woodpile  structure.  This  technique  has  several  disadvantages.  The  first  is 
the  structure  is  not  completely  metallic,  but  a  coated  structure  only.  For  our  purposes,  this  should 
not  be  a  problem  as  long  as  the  thickness  of  the  coating  is  larger  than  the  skin  depth  of  the 
incident  light.  The  second  major  problem  is  that  electroless  plating  is  not  selective.  That  means 
that  the  electroless  technique  would  plate  both  the  structure  and  the  substrate  at  the  same  time. 
To  perform  transmission  measurements  for  characterization  purposes,  it  is  important  to 
selectively  coat  the  structure  and  not  the  entire  substrate. 

Selectivity  is  a  major  problem  for  electroless  plating.  Many  techniques  have  been 
developed  to  overcome  this  problem,  mostly  by  modifying  the  substrate  surface  properties  to 
achieve  selectivity.  However,  these  techniques67  require  complicated  chemical  modification 
steps  as  well  as  long  treatment  times.  Radke  et  al’s  overcame  this  problem  using  a  rod  to  transfer 
the  silver  coated  structure  from  one  substrate  to  a  clean  substrate68.  The  transfer  part  is  the  most 
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challenging  as  structures  can  easily  break  and  degrade  in  quality  when  not  handled  properly.  In 
this  work,  we  employed  a  novel  and  simple  approach  using  radio  frequency  plasma  pretreatment 
of  the  sample  prior  to  plating  to  achieve  selectivity69. 

Electroless  plating  is  the  deposition  of  metals  on  a  catalytic  surface  which  takes  place  in 
an  aqueous  solution  without  any  external  source  of  current.  The  metal  is  reduced  from  its  ionic 
form  by  a  reducing  agent  which  is  added  into  the  solution,  before  it  is  deposited  on  the  sample 
surface.  The  quality  of  the  film  is  mainly  determined  by  the  adhesion  between  the  surface  and 
metal  nanoparticles.  In  this  work,  we  employ  the  most  widely  used  electroless  technique  based 
on  Tollens  reaction  .  Here  we  explain  the  plating  procedure  in  detail. 


1)  Prepare  water  bath:  Water  bath  is  heated  to  a  temperature  of  50°C  using  a  hot  plate.  This 
speeds  up  the  plating  process. 

2 )  Prepare  chemicals: 

a)  Silver  nitrate  solution  (AgNCE):  The  O.lmol/L  silver  nitrate  solution  is  1:1  diluted  with 
millipore  water  to  get  concentration  0.05mol/L. 

b)  Aqueous  ammonia  (NH4OH) 

c)  Reducer  solution:  Glucose  (C6Hi206) 

3)  Sample  preparation:  Surface  treatment  for  adhesion  improvement 

•  The  sample  is  immersed  in  an  acidic  solution  of  Tin  (II)  chloride  (SnCl2): 

SnCl2  :  HC1 :  CH3COOH  (2:1:3)  diluted  1:1  in  distilled  water. 

•  Immersed  for  ~  10  mins  and  thoroughly  rinsed  in  distilled  running  water. 

•  Treat  the  sample  with  RF  plasma:  dose  of  800J 

4)  Plating 


•  Add  the  diluted  silver  nitrate  solution  to  a  beaker. 

•  Add  ammonia  solution  drop  by  drop  with  magnetic  stirrer  set  at  150  rpm.  As 
ammonia  was  added,  silver  oxide  (Ag20)  precipitate  was  first  produced  and  then 
disappeared  as  more  ammonia  was  introduced  to  the  solution.  Once  the  precipitate 
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disappeared  completely,  stop  adding  ammonia.  This  is  due  to  the  formation  of 
Ammonical  silver  nitrate  complex  as  shown  below 

2AgN03  +  2NH4OH  ->  Ag20  +  2NH4N03  +  H20 

Ag20  +  4NH4OH  ->  2[Ag(NH3)2]OH  +  3H20 

•  Place  the  surface  treated  sample  onto  the  beaker  with  the  stirrer  and  add  1  ml  of 
glucose  solution.  Glucose  acts  as  the  reducing  agent  to  give  pure  silver  as  can  be 
observed  as  the  sample  gets  coated 

2[Ag(NH3)2]OH  +  CHO  ^  2Ag  +  4NH3  +  COOH  +  H20 

•  Remove  the  sample  every  5  mins  to  rinse  it  thoroughly  with  high  pressure  DI  water 
and  repeat  the  step  before  the  solution  gets  too  muddy. 

As  mentioned  earlier,  surface  roughness  is  critical  to  plasmonic  structures  and  can  affect 
the  optical  quality  of  the  structures.  However,  the  surface  roughness  of  the  third  metal  layer  is 
not  as  critical  for  our  application  as  it  is  the  metal/dielectric  interface  where  surface  plasmons 
exist  and  propagate.  The  quality  of  the  adhesion  between  the  surface  and  the  metal  coating  is 
what  we  care  the  most.  Hence,  the  plating  parameters  were  optimized  experimentally  based  on 
the  quality  of  the  metal/dielectric  interface.  Figure  2.11  shows  a  very  good  adhesion. 


Figure  2.11  -  SEM  image  of  silver  coating  on  amorphous  silicon  via  electroless  plating 
(a)  top  view  and  (b)  cross-section  view.  Image  on  the  right  also  shows  the  high  uniformity  of 
the  deposition. 
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Our  goal  is  to  use  electroless  plating  on  a  complex  3D  RGWN  structure.  To  demonstrate 
the  3D  coating  capability,  we  fabricated  silicon  inverse  opal  structure  and  coated  the  structure 
with  silver  via  electroless  plating  (Figure  2.12).  We  can  see  that  silver  is  coated  conformally  and 
uniformly  throughout  the  structure. 


Figure  2.12  -  SEM  image  of  silver  coating  on  silicon  inverse  opals  via  electroless  plating 
(a)  top  view  and  (b)  magnified  view. 


2.4  Final  structure 

After  optimizing  all  the  experimental  steps,  we  fabricated  the  3D  RGWN  structure  based 
on  a  3D  woodpile  scaffold.  In  addition,  various  factors  were  considered  while  fabricating  the 
structure  such  as  the  propagation  length  of  surface  plasmons  in  our  MIM  structure,  the 
characterization  scheme,  the  beam  spot  size  used  for  characterization... etc. 

The  schematic  of  the  final  structure  is  shown  below  (Figure  2.13).  The  structure  consists 
of  a  6-layer  3D  woodpile  with  overall  lateral  dimension  of  250  pm  x  250  pm.  The  rods  are  2.5 
pm  in  diameter  with  a  center-to-center  rod  spacing  d  of  -5.5  pm. 
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Figure  2.13  -  Schematic  of  the  final  3D  RGWN  structure. 


As  explained  later  in  Section  3.4,  we  plan  to  analyze  the  transmission  measurements  to 
study  the  optical  dispersion  within  the  3D  RGWN  structure.  In  other  words,  we  only  want  to 
detect  the  signal  that  goes  through  the  structure.  Hence,  we  infill  the  structure  with  silver  via 
electroplating  as  the  final  structure  is  a  mesh-like  structure  (Figure  2.14b).  This  eliminates  any 
background  signal  which  can  arise  from  multiple  scattering  within  the  woodpile  structure. 


Figure  2.14  -  SEM  image  of  the  final  3D-woodpile  RGWN  structure  (a)  top  view  and 
(b)  top  view  after  complete  infilling  of  the  structure  with  silver  via  electrodeposition. 
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CHAPTER  3 


OPTICAL  CHARACTERIZATION  OF  3D  RGWNs 

3.1  Introduction 

Optical  spectroscopy  is  one  of  the  most  common  techniques  to  characterize  dispersive 
photonic  materials.  Transmittance  and  reflectance  measurements  can  reveal  how  light  propagates 
and  is  manipulated  through  the  structure.  Within  the  set  of  techniques  available  for  optical 
spectroscopy,  angle  resolved  spectroscopy  is  the  most  common  and  straightforward  technique  to 
study  the  optical  dispersion  within  photonic  materials.  It  has  been  used  in  the  past  to  map  out  the 
dispersion  relations  of  2D  and  3D  photonic  crystals.  This  technique  is  ideal  to  characterize 
RGWNs  that  exhibit  wave  dispersion  and  photonic  band  gaps  due  to  interference  effects  within 
the  local  network. 

In  this  chapter,  we  introduce  realization  and  in-depth  description  of  the  main 
experimental  technique  used  to  characterize  3D  RGWNs.  Such  technique  is  based  on  direct 
inspection  of  the  spatial  optical  Fourier  transform  of  the  response  of  the  dispersive  medium  (the 
sample).  In  short,  we  will  use  angular  resolved  transmission  data  to  study  the  optical  dispersion 
within  these  structures  by  mapping  out  its  photonic  band  diagram.  We  also  present  some 
preliminary  finite  difference  time  domain  (FDTD)  simulation  results  to  validate  the  feasibility  of 
our  experimental  technique  based  on  the  theoretical  structure  constructed  from  optical 
measurements  and  scanning  electron  microscopy  (SEM). 
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3.2  Optical  response  of  the  silver  coating  (optimized) 


As  discussed  in  Section  1.3.3,  we  fabricated  a  ID  silver  line  grating  to  study  the 
plasmonic  response  of  our  final  silver  coating  since  SPPs  do  not  directly  couple  to  an  optical 
field.  We  first  printed  ID  line  gratings  with  the  silver  ink  using  direct  writing  onto  two  samples. 
One  of  the  samples  was  then  electroplated  with  silver,  electropolished  and  annealed  for 
optimization  as  discussed  in  Section  2.3.  The  grating  (Figure  3.1a)  was  coated  with  -130  nm  of 
silver  coating  which  was  later  confirmed  by  cross-sectional  SEM  image. 


Figure  3.1  -  (a)  Top  view  SEM  image  of  ID  line  grating  coated  with  silver  (optimized), 
(b)  Reflectance  spectra  of  the  plasmonic  grating  after  silver  coating  compared  to  as 
printed  Ag  nanoparticle  ink. 


Figure  3.1b  shows  a  significant  increase  in  reflectance  and  a  much  narrow  and  sharper 
resonance  peak  after  silver  coating  compared  to  the  as  printed  silver  ink.  This  enhancement  in 
the  plasmon  response  of  the  grating  is  an  indication  of  the  improved  optical  quality  of  the 
structure  after  silver  coating  in  comparison  to  the  as  printed  structure.  The  resonance  peak  shifts 
very  little  to  the  right  which  is  expected  due  to  an  increase  in  the  feature  size  of  the  grating  after 
the  coating. 
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3.3  Experimental  realization 


As  mentioned  earlier  in  Section  2.2,  the  photonic  band  structure  of  RGWNs  is  only  due 
to  dispersion  resulting  from  the  network  topology.  Hence,  we  want  to  verify  that  we  get  different 
dispersion  by  varying  the  local  network  topology  such  as  the  waveguide  length  (L)  or  the  sub¬ 
wavelength  MIM  gap  (d).  FDTD  simulations  verify  that  we  should  get  different  dispersion  with 
different  network  topology  for  such  structures  as  shown  in  Figure  3.2. 


r  xr  x 

Figure  3.2  -  Optical  DOS  of  an  infinite  3D  network  spaced  periodically  with  cubic 
periodic  unit  cell  (a)  L  =  2  pm  and  (b)  L  =  2.2  pm.  Image  reproduced  from  Reference  44. 

Figure  3.2  plots  the  optical  density  of  states  (DOS)  for  different  wave  vectors  over 
frequencies  in  the  near  infrared  range,  where  the  material  dispersion  is  small.  We  observe  a 
photonic  band  structure  which  is  only  due  to  dispersion  resulting  from  the  network  topology.  The 
k-vectors  are  extracted  from  the  frequency  spectrum  at  different  incident  angle. 

3.3.1  FIB  to  open  the  structure 

After  fabrication,  the  3D  RGWN  is  a  closed  structure.  We  open  the  structure  to  couple 

light  into  the  MIM  waveguide  structure  to  study  the  wave  dynamics  within  the  structure.  We 
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achieve  this  using  Focused  ion  beam  (FIB).  The  Focused  Ion  Beam  (FIB)  system  uses  a  Ga  +ion 
beam  to  raster  over  the  surface  of  a  sample  in  a  similar  way  as  the  electron  beam  in  a  scanning 
electron  microscope  (SEM).  The  ion  beam  allows  the  milling  of  the  sample  at  well  localized 
sites.  This  technique  is  most  widely  used  in  the  semiconductor  industry,  materials  science  to 
prepare  TEM  samples,  view  cross-sectional  images  of  a  structure  or  make  modifications  to  the 
structures. 

The  initial  idea  was  to  open  to  structure  using  FIB  cut  at  different  angles  (Figure  3.3a). 
We  then  characterize  the  samples  using  Fourier  transform  infrared  spectroscopy  (FTIR)  at 
different  FIB  cut  angles.  Amorphous  silicon  is  interconnected  throughout  the  structure  and  the 
silicon  substrate  is  transparent  in  IR.  Hence,  we  should  expect  light  output  as  we  couple  light 
into  the  structure  as  shown  in  Figure  3.3b.  This  should  allow  access  to  off-normal  incidence  k- 
vectors  to  map  out  the  dispersion  (band)  diagram  over  a  broad  range  of  k-vectors.  However,  the 
sample  is  tilted  during  the  FIB  cut  and  the  high  aspect  ratio  of  the  cut  is  very  challenging 
especially  for  angled  cuts.  Moreover,  if  the  cuts  are  not  well  defined  then  we  cannot  access  the  k- 
vectors  of  the  structure. 


Figure  3.3  -  Schematic  showing  (a)  FIB  cut  at  different  angles  to  open  the  structure, 
(b)  the  interconnected  silicon  throughout  the  structure. 
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3.3.2  Angle-resolved  Fourier  Spectroscopy 


The  usual  and  simplest  technique  for  angle  resolved  spectroscopy  consists  of  the  standard 
normal  incidence  FTIR  setup  with  two  additional  branches  mounted  on  a  goniometer  which  are 
able  to  rotate  in  a  controlled  way  around  a  sample  .  An  example  of  this  kind  of  system  is  shown 
in  Figure  3.4a.  However,  this  was  a  custom  built  setup  and  is  very  challenging  to  assemble  and 
calibrate.  Moreover,  we  only  need  transmittance  from  our  structure. 


Figure  3.4  -  Schematic  drawing  of  (a)  the  setup  for  angle -resolved  transmission  and  reflection 
optical  measurements:  A,  FTIR  spectrometer;  B,  achromatic  optical  focusing  system;  C, 
sample  and  detector  rotation  stages.  Image  reproduced  from  Reference  73.  (b)  angle-resolved 
spectroscopy  with  only  the  sample  rotated  at  defined  angles  to  access  the  off-normal  incidence 
k- vectors. 


Instead  of  doing  the  angled  cuts  via  FIB  or  use  a  complicated  setup  as  shown  above,  a  flat 
(0°)  cut  was  made  and  the  sample  was  rotated  at  defined  angles  to  access  the  off-normal 
incidence  k-vectors  (Figure  3.4b).  Note  that  this  setup  cannot  be  used  for  angle -resolved 
reflection  measurements.  The  idea  here  is  to  use  a  pinhole  setup  that  provides  a  narrow  beam 
configuration  for  the  measurement  of  photonic  band  gaps  (PBGs)  and  this  technique  has  been 
used  in  the  past  to  plot  photonic  band  structure  of  a  complex  3D  structure74.  The  drawback  of 
this  technique  is  that  as  the  plane  waves  hit  different  waveguide  at  different  times,  we  lose 
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coherence  (phase  information)  when  light  excites  the  sample  at  different  angle.  This  is  true  for 
free  space  but  the  coherence  is  maintained  at  the  point  of  incoupling  for  the  coaxial  configuration 
as  shown  in  Figure  3.5.  The  red  lines  are  the  monitors  where  the  data  was  collected  and 


analyzed.  The  plots  show  the  real  part  of  the  longitudinal  component  of  the  H-field. 
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Figure  3.5  -  Top:  Schematic  diagram  showing  the  FDTD  simulation  setup  used  to 
analyze  the  coherence  simulation  for  (a)  freespace  and  (b)  the  coaxial  MIM  geometry. 
Bottom:  Real  part  of  the  longitudinal  H-field  as  a  function  of  the  distance  light  travels 
within  a  material  (z)  for  (a)  freespace  and  (b)  for  the  coaxial  MIM  geometry. 


For  freespace,  we  can  see  that  a  1.5  pm  pulse  hits  monitor2  first  and  has  a  phase  advance 
with  respect  to  monitor  1,  thus  defining  a  k- vector.  Note  that  for  a  continuous  wave  (CW)  this 
would  be  the  other  way  around,  but  the  important  thing  here  is  that  there  is  a  phase  difference. 
Similar  dynamics  are  observed  for  a  1.5  pm  pulse  hitting  a  waveguide  array.  Thus,  the 
incoupling  does  not  cause  decoherence  of  the  incoming  phase,  allowing  for  the  structure  to  have 
a  k- vector. 
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Minor  modifications  were  made  to  a  Hyperion  2000  IR  microscope  (Bruker  Optics) 
which  was  equipped  with  Cassegrain  objectives  (15x/0.4)  and  a  liquid  nitrogen  cooled  HgCdTe 
detector  that  was  available  at  our  lab.  The  Cassegrain  objective  of  the  microscope  provides  a 
hollow  light  cone.  First,  the  numerical  aperture  (opening  angle)  of  the  15x  Cassegrain  objective 
was  reduced  which  acts  as  the  illuminating  condenser  and  sits  below  the  sample.  In  the  standard 
configuration,  the  objective  illuminates  the  sample  with  a  hollow  cone  of  wave  vectors  varying 
from  10°  to  24°  (Figure  3.6a).  Circular  disks  out  of  0.2  mm  sheet  steel  with  the  same  diameter  as 


the  top  of  the  Cassegrain  objective  were  cut  with  a  circular  hole  to  serve  as  an  aperture,  reducing 
the  opening  angle  of  the  illuminating  light  cone  to  e.g.  5°  or  10°  (Figure  3.6b). 

Standard  configuration 

a 


Figure  3.6  -  Schematic  showing  (a)  the  Standard  configuration  and  (b)  the  Narrow 
beam  configuration  for  the  Cassegrain  objective.  Image  reproduced  from  the  Newport 
website. 


For  proper  placement  of  the  holes,  we  obtained  a  technical  drawing  of  our  objective  from 
the  Newport  website.  During  the  experiment,  one  of  the  sheet  steel  disks  is  placed  on  top  of  the 
objective  as  shown  in  Figure  3.7a.  Since  the  disk  had  the  same  diameter  as  the  objective, 
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centering  it  was  easy.  As  you  can  also  see  from  Figure  3.6,  the  narrow  beam  now  hits  the  sample 
at  an  angle.  To  compensate  for  that  and  to  vary  the  angle  of  incidence  in  a  controlled  way,  a 
small  goniometer  stage  from  Thorlabs  (GN05)  and  attached  a  small  self-built  extension  arm  to  it 
for  mounting  the  sample  (Figure  3.7b).  The  sample  wasn't  actually  located  at  the  virtual  fix- 
point  of  rotation  of  the  goniometer,  so  changing  the  angle  of  incidence  would  also  move  the 
structure  laterally  relative  to  the  beam.  This  made  measurements  slightly  awkward  but  we 
couldn’t  find  a  way  around  it  in  the  limited  vertical  space  that  was  available  between  the 
illuminating  and  collecting  Cassegrain  objectives. 


Figure  3.7  -  (a)  Pinhole  setup  with  0.2  mm  sheet  steel  with  the  same  diameter  as  the 
Cassegrain  objective  on  top  of  the  objective,  (b)  goniometer  stage. 
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CHAPTER  4 


CONCLUSION  AND  FUTURE  OUTLOOK 

4.1  Conclusion 

The  search  for  new  photonic  materials  has  defined  important  directions  for  controlling 
optical  dispersion.  RGWNs  offer  a  different  approach  for  designing  dispersive  photonic 
materials.  Resonant  guided  wave  networks  are  a  new  class  of  artificial  photonic  material,  distinct 
from  photonic  crystals  and  metamaterials,  in  which  localized  waves  resonate  in  closed  paths 
throughout  a  network  of  isolated  waveguides  connected  by  power  splitting  junctions,  enabling 
the  wave  dispersion  to  be  dependent  on  the  waveguide  network  layout.  The  coherent  wave 
propagation  through  the  network  is  dependent  only  by  the  waveguide  length  and  the  phase  shift 
upon  power  splitting  and  thus  opens  up  new  directions  for  the  design  of  artificial  optical 
materials  and  devices. 

In  this  thesis,  we  first  introduce  the  underlying  physics  of  the  RGWNs  building  blocks. 
The  wave  dynamics  is  investigated  using  a  FDTD  model.  The  fabrication  of  3D  MIM  RGWN 
structures  has  been  developed.  FDTD  simulation  was  used  to  incorporate  fabrication  limitation 
to  the  model  while  ensuring  coherent  wave  propagation  and  resonant  behavior  within  these 
structures.  The  fabrication  process  is  based  on  direct  writing  of  a  scaffold  and  subsequent 
thermal  annealing,  selective  silver  plating  and  CVD  to  form  novel  3D  metallic  structures.  We 
fabricated  the  3D  RGWN  structure  based  on  optimized  experimental  procedure.  A  set  of 
experiments  to  optimize  the  plating  parameters  in  order  to  achieve  the  smallest  surface  roughness 
with  high  degree  of  control  over  the  thickness  have  also  been  performed.  We  achieved  high 
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quality  films  with  rms  roughness  ~5nm.  A  rougher  surface  leads  to  more  energy  loss  due  to 
scattering  so  it  is  not  desirable. 

We  measured  the  thickness  (~130nm,  skin  depth  ~40nm)  and  resistivity  (4.83  x  10~6  Q  • 
cm)  of  the  silver  coating  under  optimal  conditions.  FTIR  characterization  was  performed  to 
obtain  the  transmittance  spectra  of  silver  coated  ID  line  gratings.  The  results  show  good 
plasmonic  response,  verifying  the  high  quality  of  the  coating  and  feasibility  of  the  technique. 

We  then  introduce  realization  and  in-depth  description  of  the  main  experimental 
technique  used  to  characterize  3D  RGWNs.  Angle-dependent  Fourier  spectroscopy  is  used  to 
characterize  the  structure.  We  also  present  some  preliminary  finite  difference  time  domain 
(FDTD)  simulation  results  to  validate  the  feasibility  of  our  experimental  technique  based  on  the 
theoretical  structure  constructed  from  optical  measurements  and  scanning  electron  microscopy 
(SEM).  Minor  modifications  were  made  to  a  standard  Hyperion  2000  IR  microscope  (Bruker 
Optics)  to  incorporate  a  pinhole  setup  that  provides  a  narrow  beam  configuration  for  the 
measurement  of  photonic  band  gaps  (PBGs). 

In  addition,  incorporating  more  components  into  the  RGWN  could  allow  for  more 
complex  devices  or  alternatively  for  devices  with  enhanced  performance  allowing  for  the 
implementation  of  complex  optical  functions  and  could  possibly  open  new  opportunities  in  the 
design  of  photonic  circuitry  devices. 
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4.2  Future  outlook 


4. 2. 1  Fabrication 

As  discussed  in  Section  2.3.1,  the  minimum  feature  size  we  can  obtain  via  direct  writing 
is  ~2.5  |am.  The  ideal  feature  size  we  desire  is  ~1  |am  for  the  RGWN  structure  to  suppress  any 
higher  order  modes.  To  circumvent  this  problem,  we  can  use  direct  laser  writing  (DLW)  to 
fabricate  the  3D  woodpile  scaffold.  DLW  is  the  only  technique  that  allows  arbitrary  three 
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dimensional  micro  and  nanostructures  to  be  fabricated  with  sub-lOOnm  resolution  and  with  a 
relatively  high  throughput.  The  technique  utilizes  nonlinear  optical  absorption  of  a  tightly 
focused  pulsed  laser  beam  to  selectively  crosslink  a  small  volume  inside  a  photosensitive 
polymer  (Figure  4.1b).  By  moving  the  focal  point  relative  to  the  sample  in  all  three  dimensions, 
arbitrary  three  dimensional  nanostructures  can  be  fabricated  that  are  not  possible  through  other 
lithographic  technologies.  Due  to  its  high  resolution,  high  level  of  control,  low  proximity  effect, 
good  adhesion  and  easy  handling,  TPL  has  been  recently  applied  to  many  research  areas,  such  as 
photonics  and  plasmonics  ,  biomedical  microfluidics  and  MEMS  devices  .  A  standard  optical 
setup  for  direct  laser  writing  is  shown  in  Figure  4.1a. 


Figure  4.1  -  Schematic  of  (a)  a  standard  Direct-laser  writing  setup,  (b)  two-photon 
polymerization  through  a  photosensitive  material. 
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Once  we  have  a  3D-woodpile  polymer  scaffold,  fabricating  the  3D-MIM  based  RGWN 
structure  is  straightforward  since  we  have  already  optimized  all  the  experimental  procedures. 
Here,  we  propose  three  ways  to  fabricate  the  structure  and  the  cross-sectional  view  of  the  final 
MIM  unit  structure  is  shown  in  Figure  4.2. 

Option  1: 

•  Use  electroless  technique  to  coat  the  3D  scaffold  with  silver, 

•  Infill  the  structure  with  silver  via  electroplating, 

•  Open  the  structure  using  FIB,  and  burn  the  polymer 

•  Finally  infill  the  structure  with  Si  via  CVD 

Option  2: 

•  Use  the  polymer  structure  as  a  scaffold 

•  Use  electroless  technique  to  coat  both  the  metal  layers. 

Option  3: 

•  Use  positive  tone  photoresist  to  get  an  inverse  polymer  woodpile  structure 

•  Use  electroless  deposition  (bottom  up)  to  infill  the  voids,  thus  creating  the  3D 
woodpile  scaffold, 

•  Si  CVD  and  electroless  Ag  deposition. 


Option  1 

Option  2 

Option  3 

r  Si  via  CVD 

Ag  (Electroless) 

i  S*  S'  v'a  cvd 

Ag  (Electroless) 

fyCyr  Si Via  CVD 

_ ^  Ag  (Electroless  + 

Electroplating) 

Ag  (Electroless) 

Ag  (Electroless) 

Figure  4.2  -  Schematic  showing  the  cross-section  view  of  the  different  fabrication  scheme 
proposed  for  the  3D  MIM  based  RGWN  structure  via  direct  laser  writing.  Schematic  not  to 
scale. 
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In  addition  to  the  reduced  feature  size,  an  advantage  of  direct  laser  writing  is  that  the 
surface  morphology  of  the  metal  dielectric  interface  is  determined  by  the  surface  roughness  of 
the  polymer  which  is  atomically  smooth.  This  is  very  critical  for  device’s  performance  and 
eliminates  a  lot  of  surface  roughness  optimization  steps. 

4.2.2  Characterization 

Preliminary  characterization  on  test  structures  is  already  underway.  As  discussed  in 
Section  3.3,  we  propose  a  technique  to  map  out  the  photonic  band  diagram  of  the  3D  RGWN 
structure  and  study  its  wave  dynamics.  FDTD  simulation  results  validate  the  feasibility  of  the 
proposed  experimental  technique.  Using  the  modified  angled-resolved  setup,  we  plan  to  access 
the  angular  frequency  of  the  structure  at  different  off-normal  incidence  k-vectors  to  map  out  its 
photonic  band  gap.  This  is  achieved  through  accessing  transmission  spectra  of  the  structure  at 
defined  angles.  We  will  then  compare  the  experimental  results  with  simulation  to  verify  the 
novel  approach  to  optical  dispersion  using  these  resonant  structures. 
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PART  II 

PHONON  TRANSPORT  IN  COMPLEX  3D-STRUCTURE 
(SILICON  INVERSE  OPALS)* 

*This  work  was  done  in  collaboration  with  Jun  Ma  from  Sinha  Group  at  MechSE, 
University  of  Illinois  at  Urbana-Champaign. 
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CHAPTER  5 


COHERENT  PHONON-GRAIN  BOUNDARY  SCATTERING 
IN  SILICON  INVERSE  OPALS* 

*Content  in  this  chapter  was  previously  published  by  the  author  and  reproduced  with 
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Permission  .  Copyright  American  Chemical  Society,  2013. 

5.1  Background  and  general  concepts 

Thermoelectric  materials  are  promising  candidates  for  high  efficiency  power  generation 
based  on  the  Seebeck  effect  or  solid-state  cooling  devices  based  on  the  Peltier  effect.  The 
performance  of  a  thermoelectric  device  is  governed  by  the  figure  of  merit  ZT  =  S  aeff/  Cn 
, where  S,  aeff,  T,  and  Cit  are,  respectively,  the  Seebeck  coefficient,  effective  electrical 
conductivity,  absolute  temperature,  and  the  effective  thermal  conductivity.  Thermal  transport  in 
these  materials  is  a  key  factor  for  device  performance  since  the  thermoelectric  figure  of  merit  is 
inversely  proportional  to  thermal  conductivity  '  .  Nano-structured  materials  can  achieve  a 
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much  lower  thermal  conductivity  than  predicted  ,  mostly  through  phonon  boundary  scattering. 

01  oa 

As  a  result,  nano-structured  materials  such  as  single  crystal  silicon  nanowires  ’  or  membranes 
of  hole  arrays  ’  exhibit  improved  thermoelectric  performance  compared  to  their  bulk 
counterpart.  First,  as  the  size  of  structure  decreases,  its  surface  area  to  volume  ratio  increases, 
thereby  increasing  the  importance  of  boundaries  and  interfaces.  Additionally,  at  the  nanoscale 
the  characteristic  length  of  the  structure  approaches  the  phonon  wavelength,  and  other  interesting 
phenomena  such  as  dispersion  relation  modification  and  quantum  confinement  may  arise  and 
further  alter  the  thermal  transport.  Therefore,  understanding  thermal  transport  in  nanostructured 
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materials  is  crucial  for  engineering  high  performance  devices.  In  this  thesis,  we  investigate 
phonon  heat  conduction  in  silicon  inverse  opal  structures. 

R7 

Sub-micrometer  features  in  periodic  dielectric  structures  '  give  rise  to  unusual  optical 
properties  "  and  are  of  great  current  interest  in  designing  new  functionalities  in 
optoelectronics.  Phonons,  the  quanta  of  lattice  vibrations  in  a  dielectric,  can  similarly  interact  in 
interesting  ways  with  periodic  nanostructures  that  possess  features  below  the  phonon  mean  free 
path  and/or  comparable  to  phonon  wavelengths.  However,  phonon  transport  has  major 
differences  from  photon  transport  due  to  a  broadband  spectrum  and  dominant  incoherence 
arising  from  non-linearity,  in  addition  to  the  commonly  discussed  difference  in  dispersion. 
Recent  investigations  of  phonon  transport  in  nanostructures  focused  on  the  thermal 
conductivities  of  two-dimensional  “holey”  silicon91  93  and  core-shell  nanowire,94"96  hinting  at  the 
possibility  of  coherent  effects  in  phonon  transport.  Here  we  extend  the  experimental 
investigation  to  three-dimensional  periodic  nanostructures,  silicon  inverse  opals. 

The  periodicity  of  inverse  opals  influences  thermal  transport  in  two  ways:  at  continuum 
or  diffusive  length  scales,  heat  flow  lines  bend  as  they  flow  through  the  lattice  structure  resulting 
in  significant  increase  in  thermal  resistance.  At  microscopic  length  scales,  phonons  meandering 
through  the  lattice  scatter  at  surfaces  and  grain  boundaries,  reducing  the  intrinsic  thermal 
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conductivity  of  the  material.  Previous  theoretical  work  on  thermal  transport  in  inverse  opals 
has  considered  the  diffusive  regime  where  the  carrier  mean  free  path  is  much  smaller  than  the 
feature  size.  A  boundary  integral  method  solution  to  the  continuum  heat  diffusion  equation  for 
inverse  opal  geometries  shows  significant  deviation  from  the  Maxwell- Garnett  formula  for  the 
effective  thermal  conductivity  of  a  composite.  The  formula  yields  the  ratio  of  the  effective 
thermal  conductivity  to  the  material  thermal  conductivity  as  a  function  only  of  porosity  and 
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dimensionality  of  the  structure.  The  exact  solution  shows  that  the  geometry  of  the  lattice  and  the 
porosity  together  determine  the  effective  thermal  conductivity.  We  note  here  that  macroscale 
three-dimensional  periodic  structures,  termed  phononic  crystals  in  the  literature,99'100  attract 
interest  for  their  acoustic  properties  and  are  distinct  from  the  materials  investigated  here.  Our 
investigation  covers  features  on  the  order  of  10  nm,  approximately  five  orders  of  magnitude 
smaller  than  such  acoustic  periodic  structures  and  focuses  on  heat  flow  instead  of  acoustics. 

In  this  work,  we  show  that  the  lattice  periodicity  of  silicon  inverse  opals  mainly 
influences  continuum  heat  flow  and  that  microscopic  effects  on  phonon  conduction  occur 
through  scattering  at  the  grain  boundary.  The  scattering  of  phonons  from  grain  boundaries  is 
commonly  assumed  to  be  a  frequency-independent  process  across  all  temperatures.  The 
assumption  leads  to  a  phonon  mean  free  path  that  is  comparable  to  the  grain  size.  While  thermal 
conductivity  data  from  a  wide  range  of  experiments101  support  this  assumption,  recent 
experimental  work  shows  that  this  is  invalid  in  the  case  of  nanocrystalline  silicon  where  the 
scattering  rate  has  an  approximately  linear  dependence  in  frequency  ".We  report  measurements 
on  the  thermal  conductivity  of  silicon  inverse  opals  with  varying  feature  size  that  reveal  an  even 
stronger  dependence  of  the  scattering  rate  on  frequency.  We  infer  that  such  scattering  arises  due 
to  coherence  in  the  reflection  of  phonons  from  the  boundaries  of  the  inter-grain  region.  A  model 
of  thermal  conductivity  of  inverse  opals  that  incorporates  both  diffusive  thermal  transport  as  well 
as  microscopic  phonon  transport  incorporating  coherent  phonon  scattering  at  grain  boundaries, 
agrees  well  with  the  data.  We  start  below  with  details  of  fabrication  and  characterization,  and 
follow  with  thermal  conductivity  measurements  and  transport  modeling. 
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5.2  Fabrication  of  silicon  inverse  structures 


The  fabrication  of  silicon  inverse  opals  has  been  discussed  by  many  papers90'103.  Here  we 
present  an  overview  of  the  fabrication  procedure  (Figure  5.1)  that  was  used.  We  employed  the 
Stober  method  followed  by  several  regrowth  cycles  to  synthesize  silica  spheres  ranging  in 
diameter  from  300-640  nm.104  Prior  to  crystal  growth,  we  heat-treated  the  spheres  for  6  h  at  600 
°C  to  avoid  shrinkage  that  potentially  leads  to  cracking  of  the  opal  structure  during  post¬ 
processing  at  elevated  temperatures.105  Using  a  modified  vertical  deposition  method106  that 
employs  a  temperature  gradient,103  we  deposited  colloidal  crystal  films  on  a  single  side  polished 
Si  substrate  (Figure  5.1A).  After  a  Piranha  clean,  substrates  were  placed  at  a  20°  angle  in  a  20 
mL  scintillation  vial  (Fisher)  with  4  g  of  colloidal  dispersion  (2-3%  w/w  in  ethanol).  The  vials 
were  left  in  an  incubator  (Fisher,  Isotemp  125D)  at  -42-45  °C  overnight.  We  noticed  that  the 
appropriate  temperature  window  used  to  grow  high  quality  opals  was  small  and  occasionally 
needed  to  be  experimentally  adjusted.  Reflection  and  transmission  spectroscopy  along  with  SEM 
provided  measures  of  the  thickness  and  quality  of  the  samples.  We  only  selected  high  quality 
(closely-packed)  templates  with  -5-8  opal  layers  for  further  processing.  Figure  5.2A  shows  an 
SEM  of  a  typical  sample. 

For  such  templates,  we  deposited  a  thin  conformal  film  of  high-quality,  low-roughness 
amorphous  Si  around  the  spheres  (Figure  5.1B)  using  a  static  CVD  system  with  disilane  (SEH6, 
98%,  Gelest)  as  a  silicon  source  at  350  °C  (-50  mbarr,  heating  rate  8  °C  min'1).  We  controlled 

1  ryj 

the  deposition  time  to  make  sure  that  the  pores  in  the  opal  template  were  filled  to  a  maximum. 
Pores  in  the  underlying  opal  structure  are  filled  to  the  maximum  at  86%  once  the  silicon 
overlayer  starts  forming  during  the  fabrication  described  above.  The  shell  thickness  is 
approximately  equal  to  6%  of  the  diameter  of  silica  in  this  case.109  Scanning  electron 
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microscopy  (SEM)  provided  additional  confirmation  of  this  feature  size.  Samples  were  annealed 
in  a  tube  furnace  at  1000  °C  for  10  hours  under  forming  gas  to  recrystallize  the  amorphous 
silicon.  We  used  reactive  ion  etching  (1  min,  70  W,  gasses  SF6  and  CE  20  seem,  50  mTorr 
chamber  pressure)  to  open  a  1x2  mm  window  defined  by  a  Kapton  film  (Figure  5.1C).  This 
window  allows  the  creation  of  an  interpenetrating  network  of  oxide  within  the  opal  template  by 
exposing  the  silica  spheres  at  the  surface  when  immersed  in  an  ethanolic  solution  of  5% 
hydrofluoric  acid  for  -25-30  mins.  Selectively  etching  only  a  small  region  of  the  sample 
prevented  cracking  and  sample  lift-off.  The  opal  template  was  completely  etched  away  leaving 
the  silicon  inverse  opal  structure  (Figure  5.1D).  Figure  5.2B  shows  the  SEM  of  a  typical  sample 


after  RIE  and  HF  etching. 

A 


B 


Opal  Template 


Si  CVD 

o 


Si  Inverse  Structure  Anneal/RIE 

Figure  5.1  -  Overview  of  the  silicon  inverse  opal  structure  fabrication  scheme.  A.  3D  fee 
silica  opal  template.  B.  Silica  opal  coated  with  a  thin  layer  of  a-Si  via  CVD.  C.  Annealed 
sample  with  RIE  to  expose  silica.  D.  Silicon  inverse  structure  after  etching  the  silica 
template. 
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5.3  Shell  thickness  and  grain  size  characterization 


The  diameter  of  silica  spheres,  D  controls  two  feature  sizes  in  the  inverse  opal  that 
directly  affect  phonon  heat  conduction:  The  thickness  of  the  silicon  shell,  d,  that  controls  the 
phonon  mean  free  path  for  surface  scattering  and  the  nominal  size  of  the  grains,  /,  that  controls 
phonon  grain  boundary  scattering.  We  have  discussed  measurement  of  the  shell  thickness  above 
and  now  discuss  measurements  of  the  grain  size.  We  anticipate  grains  to  be  limited  in  size  by  the 
thickness  of  the  shell.110  Transmission  electron  microscopy  (TEM)  confirms  polycrystallinity 
and  also  reveals  atomic  scale  surface  roughness  but  can  only  image  the  grain  over  a  small  region. 


Figure  5.2  -  A.  SEM  image  of  the  opal  template.  B.  Top  view  SEM  of  the  structure  after 
RTE  opening  and  HF  etching.  C.  TEM  showing  roughness  scale  on  the  surface  of  the  silicon 
shell.  Inset  shows  selected  area  diffraction  pattern.  D.  High  resolution  TEM  showing  the 
polycrystalline  nature  of  the  sample. 


Figure  5.2C  and  D  show  TEMs  of  a  typical  sample.  We  employed  X-ray  diffraction 
(XRD)  to  measure  the  average  grain  size  across  the  entire  depth  of  the  sample  over  spot  sizes 
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~lmm  in  diameter.  A  2#-&>  scan  yields  the  average  crystallite  size  along  the  height  of  the 
inverse  opal  layer.  The  inset  of  Figure  5.3  shows  the  XRD  data  for  a  sample  fabricated  using 
420  nm  silica  spheres.  Figure  5.3  further  plots  the  average  grain  size  for  the  measured  inverse 
opals  as  a  function  of  the  diameter  of  the  silica  spheres  used  in  the  fabrication.  The  systematic 
trend  of  the  data  confirms  the  hypothesis  that  grain  size  is  essentially  limited  by  the  thickness  of 
the  silicon  shell.  As  discussed  above,  the  diameter  of  the  silica  spheres  directly  controls  the 
thickness  of  the  silicon  shell.  We  find  the  average  grain  size  to  decrease  from  36  nm  (with 
D=640  nm)  to  20  nm  (with  D=300  nm).  As  the  diameter  of  the  spheres  and  consequently  the 
thickness  of  the  silicon  shell  decreases,  the  distribution  in  grain  sizes  also  becomes  narrower. 


Figure  5.3  -  The  average  grain  size  from  XRD  measurements.  Inset:  XRD  profile 
of  a  typical  silicon  inverse  opal  sample.  We  average  the  grain  size  over  three 
orientations,  (111),  (220)  and  (311). 
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5.4  Thermal  conductivity  measurements 


5.4.1  The  3 co  method 

Any  measurement  on  thermal  conductivity  requires  two  quantities:  the  heat  flow  and  the 
temperature  gradient.  Our  samples  effectively  possess  a  thin-film  geometry  with  thicknesses  in 
the  range  2-6  pm.  The  frequency  domain  3 co  method111  can  readily  measure  the  thermal 
conductivity  of  such  samples.  A  3 co  device  usually  consists  of  a  bulk  substrate,  a  thin  film  of 
interest  and  a  heater  (thermometer)  as  shown  in  Figure  5.4.  Previous  work  ’  describes  the  3 co 
method  in  detail.  In  brief,  a  sinusoidal  current  at  a  frequency  co  through  a  metallic  heater  sets  up 
an  oscillatory  temperature  rise,  AT  at  the  harmonic  frequency  2 co.  This  results  in  a  voltage  across 
the  metal  line  oscillating  at  frequency  3 co.  Measurement  of  the  3 co  voltage  via  a  lock-in  amplifier 
provides  AT.  Comparing  the  measured  AT  with  a  model  of  multilayer  heat  diffusion114  enables 
the  extraction  of  thermal  conductivity  of  the  sample.  We  next  discuss  the  design  of  3 co 
measurements  for  our  samples. 


Top  view 

Figure  5.4  -  Schematic  of  a  3 co  device. 

Based  on  previous  measurements  of  thin  silicon115  and  polysilicon  films116,  we  expect  the 
effective  room  temperature  thermal  conductivity  of  the  films  to  be  rather  small  at  ~1  W/mK  on 


Metal  heater 

Film 

Silicon  substrate 

Cross  section  view 
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account  of  the  large  porosities.  Assuming  a  thermal  conductivity  of  1  W/mK  and  a  bulk-like 
volumetric  heat  capacity  scaled  appropriately  by  the  porosity,  the  expected  thermal  diffusivity,  a 

is  ~  3x10  2cm2/s  at  room  temperature.  For  the  typical  frequency  range  of  10-1000  Hz  used  in 
3 co  measurements,  we  estimated  the  thermal  penetration  depth  to  be  -30-300  /./ m  at  room 
temperature.  This  is  much  larger  than  the  thickness  of  the  inverse  opal  layer  and  allows  the  3 co 
measurement  to  probe  the  entire  thickness  of  the  thin  film  with  the  substrate  underneath  acting  as 
the  reference.  Heat  flow  across  the  sample  during  the  measurement  is  essentially  one¬ 
dimensional  when  the  width  of  the  heater  is  chosen  to  be  greater  than  the  thickness  of  the 
sample  .  This  further  simplifies  thermal  conductivity  extraction. 

The  preparation  of  samples  for  the  3 co  measurements  proceeded  as  follows.  We  employed 
plasma-enhanced  chemical  vapor  deposition  (PECVD)  to  deposit  a  100  nm  thick  silicon  dioxide 
layer  on  the  sample,  and  patterned  a  300  nm  thick,  50  /r m  wide  Au  heater  line  through  a  shadow 
mask.  The  heater  is  thick  enough  to  maintain  physical  contact  with  the  wavy  substrate 
underneath.  The  large  width  compared  to  the  sample  thickness  ensures  one-dimensional  heat 
flow  across  the  sample,  as  discussed  above.  Post-measurement  cross-sectional  SEMs  confirmed 
continuous  physical  contact  across  all  interfaces. 

5.4.2  Data  and  discussion 

All  layers  and  interfaces  present  in  the  3 co  sample  described  above  contribute  to  AT  but 
the  thermal  impedance  posed  by  the  inverse  opal  layer  dominates  over  the  rest.  The  thermal 
diffusivity  of  the  bulk  Si  substrate  is  typically  more  than  two  orders  of  magnitude  larger  than  that 
of  the  inverse  opal  layer.  At  300  K,  we  expect  the  inverse  opal  layer  to  dominate  and  contribute 
more  than  80%  of  the  temperature  rise  at  the  heater.  Based  on  previously  measured  values  ,  we 
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estimate  impedances  from  interfaces  and  the  oxide  to  be  small.  At  300  K,  for  example,  these 
contribute  <3%  to  the  total  temperature  rise.  Calibration  for  the  measurement  involves 
measurement  of  the  electrical  resistance  of  the  heater  line  at  low  currents  (-250  nA).  We  use  the 
Bloch-Griineisen  formula  to  fit  the  electrical  resistance  of  the  heater  line  in  the  temperature 
range  of  15-400  K.  This  serves  as  the  calibration  for  subsequent  electrical  resistance  based 
thermometry.  Differentiating  the  fit  yields  the  temperature  coefficient  of  resistance,  dR  /  dT  of 
the  heater. 

During  the  measurements,  we  observed  a  linear  relation  between  AT  and  the  logarithm  of 
heater  frequency  for  all  temperatures,  confirming  an  essentially  one -dimensional  heat  flow 
across  the  inverse  opal  films.  We  checked  for  errors  in  the  one-dimensional  approximation  by 
fitting  the  temperature  rise  with  a  two-dimensional  multilayer  heat  diffusion  model  but  found  the 
discrepancy  in  extracted  thermal  conductivities  to  be  <0.5%.  We  note  that  during  the  course  of  a 
complete  temperature  sweep  (30-400  K),  AT  was  always  smaller  than  2  K  for  all  samples  except 
for  the  smallest  sample  at  30  K.  In  the  latter  case,  AT  was  -4  K.  Since  we  assume  the  base 
temperature  of  the  sample  as  the  temperature  of  the  cryostat,  we  expect  the  error  in  this 
assumption  to  be  less  than  2  K  except  as  noted  above.  We  could  not  obtain  reliable  data  below 
30  K  since  the  resistivity  of  Au  loses  sensitivity  to  temperature. 

Figure  5.5  shows  the  measured  effective  thermal  conductivities  (right  axis)  of  four 
inverse  opal  samples119.  The  thickness  of  the  silicon  shell  and  the  average  grain  size  in  each 
sample  are  indicated.  All  samples  possess  a  relatively  low  thermal  conductivity  (-1  W/mK  at 
room  temperature)  owing  to  the  large  porosity  and  bending  of  heat  flow  lines  inside  the  structure 
as  discussed  previously.  Using  numerical  solutions  obtained  in  Ref.  ,  we  convert  the  effective 
thermal  conductivities  to  the  material  thermal  conductivities.  The  latter  are  shown  on  the  left 
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axis.  The  thermal  conductivity  of  previously  measured  undoped  LPCVD  poly  silicon  films  with 
~l-4  ju m  thickness101  ranges  between  -15-55  W/mK  at  300  K  for  grain  sizes  in  the  range  -190- 
550  nm.  Thermal  conductivity  modeling  of  this  data  set  revealed  that  thermal  conductivity  scales 
approximately  proportional  to  the  grain  size.  The  material  thermal  conductivities  of  silicon 
inverse  opals  are  lower  than  those  of  the  previous  LPCVD  thin-films  since  the  grain  size  in  our 
samples  are  nearly  an  order  of  magnitude  smaller  than  previous  samples.  However,  in  contrast  to 
LPCVD  thin  films,  the  material  thermal  conductivities  of  inverse  opals  at  300  K  do  not  scale  in 
direct  proportion  to  the  grain  size  when  compared  with  the  previously  measured  values  or 
between  our  samples.  A  direct  scaling  of  thermal  conductivity  with  grain  size  would  yield 
material  thermal  conductivities  of  inverse  opals  in  the  range  2.8-5. 1  W/mK  at  300K.  The 
measured  values  are  approximately  50%  larger.  Thus,  a  grain  size  dependent  mean  free  path  does 
not  explain  the  thermal  conductivity  of  inverse  opals.  We  now  turn  to  the  low  temperature  data 
for  more  insight. 


Figure  5.5  -  Measured  effective  (right  axis)  and  material  (left  axis)  thermal  conductivities 

of  silicon  inverse  opals.  The  fitted  curves  are  calculated  from  the  modified  Callaway's 

2 

model  with  -  (0  dependence  in  the  phonon  grain  boundary  scattering  rate.  The  error  bars 
on  each  curve  indicate  the  percentage  uncertainty  of  that  sample. 
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The  variation  of  thermal  conductivity  with  temperature  at  low  temperatures  (prior  to  the 
strong  onset  of  Umklapp  scattering)  is  indicative  of  the  frequency  dependence  of  the  dominant 
phonon  scattering  mechanism.  Applying  results  from  the  kinetic  theory,  the  lattice  thermal 

1 90 

conductivity  of  a  three-dimensional  phonon  gas  is 

k  =  \LCv A  (1) 

j,a> 

where  j  represents  different  phonon  modes.  The  volumetric  heat  capacity,  C,  the  phonon  group 
velocity,  v,  and  the  phonon  mean  free  path,  A  are  frequency  dependent  quantities  appearing  in 
the  conductivity  summation.  At  low  temperatures,  if  the  surface  or  alternately  the  grain  boundary 
dominates  phonon  scattering,  the  mean  free  path  typically  remains  independent  of  the  phonon 
frequency.  Thus,  the  temperature  dependence  of  thermal  conductivity  follows  the  well-known  ~ 
T3  dependence  of  the  heat  capacity.  The  thermal  conductivity  of  LPCVD  polysilicon  films 
follows  such  temperature  dependence  at  low  temperatures  indicating  a  frequency-independent 
mean  free  path  arising  from  grain  boundary  scattering. 

The  conductivity  data  in  Figure  5.5  however,  varies  as  k~T18  at  low  temperatures.  We 
expect  the  material  in  the  sample  to  possess  a  heat  capacity  similar  to  the  bulk  since  lattice 
dynamics  calculations  clearly  show  phonon  dispersion  and  density  of  states  to  remain  unaltered 
at  these  dimensions  “  .  In  silicon  inverse  opals,  surface  scattering  is  equally  important  as  grain 
boundary  scattering  since  the  feature  sizes  are  similar.  However,  a  frequency  dependence  arising 
from  surface  scattering  is  unlikely  since  TEM  images  show  atomic  scale  roughness  at  the 
surfaces.  In  this  case,  surface  scattering  should  be  consistent  with  the  Casimir  model  ’  and 
independent  of  frequency.  Thus,  the  weaker  than  ~  T  '  temperature  trend  of  the  data  is  puzzling  at 
first  glance.  Using  detailed  phonon  scattering  theory  below,  we  show  that  this  trend  plausibly 
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originates  from  coherence  of  the  reflected  phonon  in  the  inter-grain  region.  This  yields  a 
frequency  dependence  in  the  grain  boundary  scattering  rate  that  results  in  a  weaker  than  ~T 3 
dependence  in  the  thermal  conductivity  data. 


5.4.3  Phonon-grain  boundary  scattering  theory 


In  equating  the  phonon  mean  free  path  with  the  grain  size,  there  is  an  implicit  assumption 
that  phonons  are  randomly  scattered  at  all  grain  boundaries.  This  is  not  strictly  valid  at  low 
temperatures.  Phonons  scatter  near  a  grain  boundary  due  to  change  in  orientation  of  crystal 
planes  as  well  as  disorder  in  the  inter-grain  region.  In  a  perturbative  approach,  both  induce  a 
perturbation  in  the  wave  velocity,  Av.  The  change  in  wave  velocity  leads  to  back-scattering  from 
mode  q  to  q’  whose  rate  is 


1 

r,(q) 


V 

(2^7 


p-V22§v;12|M(q,q')f^' 

j' 


(2) 


where  V  is  the  volume  of  a  crystal  containing  one  boundary,  p  is  the  density,  vg  is  the  phonon 

velocity,  dS  is  the  surface  element,  j  is  polarization.  The  perturbation  element  )  is 

M (q,  q ')  =  V1  Ipco1  J drSv(r)e,Qr  fe  -  s')  (3) 

124 

where  £,  s’  are  unit  vectors  along  the  polarization  directions,  and  Q=q’-q.  Klemens  showed 
that  randomized  scattering  at  a  tilt  boundary  gives  rise  to  a  frequency  independent  scattering  rate 

7-  =  |(v//)(  Av/v)2  (4) 
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where  /  is  average  grain  size,  Av  is  the  change  is  velocity  at  the  tilt  boundary.  The  magnitude  of 
the  rate  is  typically  small  since  Av/v  ~  o(0.1).  Similar  randomized  scattering  in  the  inter-grain 
region  with  greater  disorder  leads  to  much  stronger  scattering  since  Av/v  ~  1  in  this  case. 
However,  for  scattering  to  be  incoherent,  the  two  boundaries  of  the  inter-grain  region  should  be 
far  apart  compared  to  the  phonon  wavelength.  This  suggests  a  critical  frequency, cocr=v/t , 
below  which  scattering  is  increasingly  coherent.  Here,  t  is  the  thickness  of  the  inter-grain  region. 

For  frequencies  lower  than  this  critical  frequency,  the  two  boundaries  scatter  phonons 
coherently.  This  is  especially  likely  at  low  temperatures  since  t  is  typically  ~A°  scale  and  much 
smaller  than  phonon  wavelengths.  The  scattering  rate  in  this  case  has  a  quadratic  frequency 
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dependence  "  , 
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where  Av  is  now  the  change  in  velocity  due  to  inter- grain  region.  We  show  below  that  such 
frequency  dependence  explains  the  weaker  than  ~T}  dependence  in  the  thermal  conductivity  of 
inverse  opals. 

We  also  point  out  recent  work  on  nanocrystalline  Si  that  suggests  linear  frequency 
dependence  in  the  grain  boundary  scattering  rate.  The  expression  for  the  scattering  rate  in  this 
case  is  proportional  to  the  phonon  transmission  coefficient  of  the  grain  boundary.  Atomistic 
simulations  reveal  that  the  transmission  coefficient  itself  can  have  strong  frequency 
dependence,  in  which  case  the  frequency  dependence  of  the  scattering  rate  would  be  stronger 
than  the  linear  dependence  assumed  in  the  work.  To  proceed,  we  fit  the  thermal  conductivity  data 
using  the  ~a>°,  ~co\  -of  frequency  dependencies  respectively  in  the  grain  boundary  scattering 
rate  discussed  above. 
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5.4.4  Thermal  conductivity  modeling 


To  obtain  the  fit,  we  follow  an  approach  to  modeling  thermal  conductivity  similar  to 
those  proposed  by  Callaway  and  Holland  “  but  include  a  modification  proposed  by  Mingo  “  . 
Mingo  suggests  the  use  of  a  cut  off  frequency,  coc  smaller  than  the  Debye  frequency,  cod  in  fitting 
the  experimental  thermal  conductivity  of  silicon  nanostructures.  This  approach  accounts  for  the 
difference  in  the  phonon  dispersion  relation  between  the  Holland  model  and  more  accurate 
lattice  dynamics  calculations.  Here,  we  follow  Mingo’s  approach  and  treat  coc  as  a  fitting 
parameter  in  subsequent  modeling.  We  note  that  the  fit  is  not  overly  sensitive  to  the  particular 
value  of  coc  since  Umklapp  scattering  strongly  reduces  the  contribution  of  high  frequency 
phonons  to  heat  conduction  much  before  the  cut-off  frequency  '  . 

Assuming  a  linear  dispersion  relation  and  treating  all  polarizations  equivalently,  the 
expression  for  thermal  conductivity  is 


k  = 


2  n-kBT\ 
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where  v  =  6.4xl03m/s  is  now  the  average  phonon  propagation  speed127,  kB  is  the  Boltzmann 
constant  and  h  is  the  reduced  Planck's  constant.  The  frequency  dependent  phonon  relaxation 
time  z(co)  can  be  calculated  using  Matthiessen’s  rule,  assuming  scattering  due  to  isotopes, 
Umklapp  processes,  surfaces  and  grain  boundaries  to  be  independent.  We  use  scattering  rates 
from  the  literature126'129  for  mass  difference  (isotope)  and  Umklapp  scattering.  These  are  of  the 

form  rm  =  Aco  and  ru  =  BT ore  .  The  value  of  A  depends  on  the  isotope  concentration  and 
we  assume  the  commonly  used  value  ,  A  =  1.32xl0~45  s  .  We  follow  Ref.  in  determining 
the  values  of  B ,  C  and  eoc .  This  approach  yields  fits  to  the  high  temperature  thermal 
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conductivities  of  bulk  silicon  and  inverse  opals  using  identical  values  of  the  parameters, 
5  =  1.6x1019s/K,  C  =  152  K  and  coc  =  4.3xl013  rad/s.  We  note  that  the  cut-off  frequency 
obtained  here  compares  well  with  Mingo’s  fits  to  single  crystal  nanowires,  which  yields 
4.2  xlO13  rad/s. 

Surface  scattering  in  an  inverse  opal  is  complicated  due  to  the  complex  geometry.  A 
rigorous  but  tedious  approach  involves  solving  the  phonon  Boltzmann  equation  in  the  exact 
geometry.  However,  the  main  effect  of  the  lattice  geometry  is  the  macroscopic  bending  of  heat 
flux  lines  rather  than  microscopic  effects  on  phonon  transport,  due  to  the  relatively  large 
curvatures  involved.  Further,  grain  boundaries  likely  randomize  phonon  scattering  and  remove 
any  geometrical  effect  of  the  lattice  in  surface  scattering.  With  this  in  mind,  we  adopt  a  simpler 
approach  in  modeling  surface  scattering.  We  assume  scattering  from  the  surfaces  of  the  thin 
segments  of  the  inverse  opal  to  be  the  same  as  scattering  from  the  surface  of  a  thin  film.  Using 
Casimir’s  approach,  we  obtain  rh  1  =v/(Fd )  where  F= 4  for  a  thin-film  geometry130.  We  show 
below  that  this  approach  leads  to  excellent  fits  for  all  samples,  validating  our  simplifying 
assumptions. 

5.5  Results  and  discussion 

Figure  5.6  plots  the  scattering  rates  obtained  from  best  fits  to  the  data.  Surface  scattering 
appears  as  a  flat  line  in  the  Casimir  approximation  and  is  based  on  the  shell  thickness  as  the 
feature  size.  Since  the  grain  size  in  these  samples  is  approximately  equal  to  the  shell  thickness, 
the  surface  scattering  rate  is  approximately  the  same  as  that  for  frequency  independent  grain 
boundary  scattering.  Overall,  surface  scattering  dominates  at  low  frequencies  and  grain  boundary 


61 


scattering  dominates  at  higher  frequencies.  The  contribution  from  grain  boundary  scattering  is 
always  significant  except  for  very  low  frequencies  that  contribute  little  to  heat  conduction  across 
the  temperature  range  of  the  data.  The  previously  suggested  ~o)  model  and  the 
coherent  ~  of  models  yield  rates  larger  than  surface  scattering,  for  phonon  frequencies  exceeding 
~1  THz.  The  ~o)  rate  diverges  from  the  ~oj  rate  at  frequencies  lower  than  ~1.5  THz  and  is 
approximately  one  order  of  magnitude  smaller  at  the  lowest  frequencies.  We  note  that  Umklapp 
scattering  rates  are  more  than  an  order  of  magnitude  smaller  than  the  dominant  rates,  even  at  300 
K.  Overall,  these  calculations  suggest  that  the  contribution  of  a  frequency  dependent  grain 
boundary  scattering  alters  the  temperature  dependence  of  the  thermal  conductivity  data  from  the 

~T3  dependence  when  lower  phonon  frequencies  dominate  heat  conduction.  The  frequency 
dependence  is  insignificant  when  the  higher  frequency  part  of  the  spectrum  in  Figure  5.6 
dominates  heat  conduction,  such  as  at  room  temperatures. 


Phonon  Frequency,  (o/2tt  (THz) 

Figure  5.6  -  Frequency  dependent  phonon  scattering  rates  calculated  for  the 
inverse  opal  with  smallest  features  (d=  1 8  nm,  1=20  nm).  Rates  for  Umklapp 
scattering  are  shown  at  30  K  and  300  K.  The  deviation  between  two  grain  boundary 
scattering  models  becomes  significant  at  frequencies  below  2  THz. 
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Referring  back  to  Figure  5.5,  all  three  models  for  grain  boundary  scattering  fit  the  data 
reasonably  well  at  temperature  above  200  K.  Figure  5.5  only  shows  the  fit  obtained  using 
the  ~  a}  frequency  dependent  model  but  the  ~  oj\  ~  co  models  fit  the  data  equally  well  near  room 
temperature.  However,  the  low  temperature  behavior  is  significantly  different.  Figure  5.7  shows 
the  low  temperature  data  with  error  bars  and  includes  the  best  fits  from  the 
~a>°,  ~a>1,  ~o}  scattering  rates.  We  find  that  only  the  ~of  dependent  rate  corresponding  to 
coherent  scattering  in  the  inter-grain  region,  is  able  to  fit  the  data  reasonably  well  at  the  lowest 
temperatures.  There  is  a  systematic  trend  in  the  agreement  of  the  other  models  with  the  data. 
Predictions  from  the  frequency  independent  model  increasingly  diverge  from  the  data  below 
temperatures  of  -200  K  whereas  predictions  from  the  ~co  model  increasingly  diverge  from  the 
data  below  -140  K.  Only  the  ~  of  model  fits  the  data  down  to  30  K.  Further,  the  values  of  the 
inter-grain  thickness  obtained  from  the  best  fits  range  between  1.8-2. 2  A0,  consistent  with  values 
calculated  using  molecular  dynamics  of  bulk  silicon  for  grain  boundaries  of  different 

IOC  ni 

energies  ’  .In  fact,  the  frequency  dependence  of  the  phonon  reflection  coefficient  calculated 
from  these  simulations  follows  ~o2'3  2  5,  approximately  consistent  with  the  proposed  model.  The 
excellent  agreement  of  the  data  with  the  coherent  model  and  the  agreement  of  the  frequency 
dependency  in  the  model  with  atomistic  calculations,  strongly  suggest  that  grain  boundary 
scattering  is  indeed  coherent  in  silicon  inverse  opals  at  low  temperatures.  The  reason  behind  this 
is  not  clear  at  present.  An  obvious  hypothesis  is  that  grain  boundaries  in  silicon  inverse  opals 
possess  relatively  lower  disorder  and/or  are  thinner  compared  to  bulk  or  thin  film  polysilicon.  In 
either  case,  this  enables  coherence  in  phonon  reflections  from  the  inter-grain  region  at  low 
temperatures. 
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Temperature  (K) 

Figure  5.7  -  The  material  thermal  conductivities  of  inverse  opals  vary  as  at  low 

temperatures.  Quadratic  frequency  dependence  in  the  phonon  grain  boundary  scattering 
rate  yields  excellent  agreement  with  data.  The  solid  curves  use  ~  or  dependent  grain 
boundary  scattering  rate  while  the  dashed  and  dotted  curves  use  ~  o)  and  ~  co°  dependent 
rates.  The  ~  or  dependence  arises  due  to  coherent  scattering  in  the  inter-grain  region. 


5.6  Conclusion 

In  conclusion,  thermal  conductivity  measurements  on  silicon  inverse  opals  show  effective 
thermal  conductivities  551  W/mK  and  material  thermal  conductivities  5510  W/mK  at  room 
temperature.  The  relatively  low  thermal  conductivities  are  significant  in  photonic  applications 
where  even  relatively  small  absorption  can  cause  significant  temperature  rise.  In  exploring 
phonon  transport  in  these  structures,  we  find  that  frequency-dependent,  coherent  phonon  grain 
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boundary  scattering  explains  the  data  across  the  temperature  range  of  the  measurement  down  to 
30  K.  Compared  to  previous  measurements,  inverse  opals  provide  access  to  more  uniform  grains 
across  the  entire  sample  since  grain  growth  is  restricted  by  the  thickness  of  the  silicon  shell. 
Based  on  the  thermal  conductivity  analysis,  we  hypothesize  that  the  inter-grain  region  is  thinner 
and  possesses  lesser  disorder  than  typical  polysilicon  films.  This  is  likely  the  reason  for  the  clear 
frequency  dependence  in  our  data  at  low  temperatures  that  is  not  observed  in  previous 
measurements.  Assuming  that  inverse  opals  can  be  heavily  doped  similar  to  bulk  polysilicon, 
these  materials  become  interesting  for  thermoelectric  energy  conversion  at  high  temperatures. 

i  in 

Theoretical  calculations  “  show  ZT-0.6  at  600  K,  providing  impetus  for  future  experiments. 
This  work  provides  thermal  conductivity  data  useful  in  technological  applications  of  inverse 
opals  and  insight  into  the  physics  of  phonon  heat  conduction  in  these  structures. 
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